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Abstract 

We prove local existence of solutions for a class of suitably renormalised coupled 
SPDE-ODE systems driven by space-time white noise, where the space dimension 
is equal to 2 or 3. This class includes in particular the FitzHugh-Nagumo system 
describing the evolution of action potentials of a large population of neurons, as well 
as models with multidimensional gating variables. The proof relies on the theory of 
regularity structures recently developed by M. Hairer, which is extended to include 
situations with semigroups that are not regularising in space. We also provide explicit 
expressions for the renormalisation constants, for a large class of cubic nonlinearities. 
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1 Introduction 

Many spatially extended systems can be described by a reaction-diffusion equation coupled 
to an ordinary differential equation, of the form 

dtu = T)A x u + F(u, v ) , 

d t v = G(u,v ), (1.1) 

where u = u(t, x) and v = v(t, x ) are functions of time t ^ 0 and space x G R rf and take 
values in R m and R n respectively, A x denotes the usual Laplacian, (D E R mxm - g a p OS itfv e 
definite matrix of diffusion coefficients, and F : R m x R n —> R m and G : R m x R n —> R n 
are given sufficiently smooth functions. A well-studied situation of this kind arises in 
neuroscience, where i£R measures the position along a neuron’s axon, u(t, x) € R is the 
membrane potential at time t and position x, and the so-called gating variables v G R n 
describe the number of open ion channels of n different types. Classical examples are the 
Hodgkin-Huxley model [24] and its simplifications such as the Morris-Lecar model [32] and 
the FitzHugh-Nagumo model [15(133], However, there are many more interesting situations 
outside neuroscience that can be described by coupled systems of the form (11.111 . for 
instance in chemical kinetics [dO], in population dynamics [8], and in pattern formation [2], 
Furthermore, many other classical reaction-diffusion equations from mathematical biology, 
such as the Keller-Segel [25] or the Gray-Scott model m are of the form m if one regards 
the model parameters v not as constants via dtv = 0 but as (slowly-)changing parameters 
via dtv = £G(u, v) for some small parameter 0 < e 1 [[28] and some smooth vector field 
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G. In fact, the number of applications of systems of the form (11.111 goes far beyond the 
areas we mentioned here. 

However, in some cases, a deterministic model of the form a is not sufficient to 
describe the dynamics, and one has to add a stochastic term to capture the effect of thermal 
fluctuations and other unresolved dynamical processes affecting the system. In the case 
of neuron dynamics, see for instance the book and recent survey by Tuckwell [Ml EI| and 
the recent reviews by Bressloff 6, 7J. In particular, stochastic versions of the FitzHugh- 
Nagurno PDE have attracted considerable recent attention mmm as well as related 
models of stochastic neural fields [29|f27j. In the derivation of these models and in related 
numerical studies, one frequently uses correlated [31] as well as space-time white noise |lj 
stochastic forcing. 

Before studying properties of solutions to a stochastic partial differential equation 
(SPDE) with space-time white noise, one has to ensure that such an equation is well- 
defined, i.e., one has to attach a meaning of “solution” to a given SPDE. Depending on 
the space dimension d, this problem may in fact be extremely hard. Consider the case of 
a one-component system (m = 1, n = 0, (9 = 1) of the form 

d t u = A x u + F{u) + £ , (1.2) 

where £ denotes space-time white noise (precise mathematical definitions will be given 
below). If d = 1, it is possible to define a notion of mild solution via the Duhamel 
principle in a quite general setting (see for instance [!!!)■ In higher space dimension, the 
problem is much more difficult, owing to the fact that space-time white noise is extremely 
singular. 

Consider for instance the case of the Allen-Cahn equation given by (11.21) with F{u ) = 
au — u 3 (if a = 0 this is also known in quantum held theory as the dynamical model, 
cf [13] '). In the case d = 1, the proof of existence of a unique solution goes back to Faris 
and Jona-Lasinio m, and many quantitative properties of this solution are known (see 
for instance 0® [3]). The two-dimensional case was solved by Da Prato and Debussche 
using a particular class of Besov spaces [10] . The case d = 3, however, was only solved very 
recently by Hairer, using his theory of regularity structures |I9j. One of the difficulties 
is that a renormalisation procedure has to be used to properly define solutions. This 
is achieved by replacing space-time white noise £ by a mollified version £ £ , solving the 
resulting regularised equation, and passing to the limit of vanishing regularisation (for an 
alternative approach, see [30]). 

The theory of regularity structures provides a framework to study SPDEs with very 
singular noise, including but not limited to parabolic equations of the form (11 .2]) . The 
basic idea is to construct an abstract space in which one can define algebraic operations 
on distributions (in the sense of generalised functions), such as multiplication, composition 
with a smooth function, and convolution with a kernel. The fixed-point equation obtained 
by applying Duhamel’s principle to the SPDE with mollified noise is then lifted to the ab¬ 
stract space, solved in that space, and finally projected down to a distribution in “physical” 
space. In addition, the theory allows to incorporate the renormalisation procedure that 
is needed in most cases; see for instance [18] for an introduction. Furthermore, there are 
already several recent applications of regularity structures, e.g., to large deviation theory 
of the Allen-Cahn equation [23], to the dynamical sine-Gordon model [22]. to Wong-Zakai 
approximation of nonlinear parabolic SPDEs [21], to the parabolic Anderson model on 
bounded [19] as well as unbounded m domains, and to the KPZ equation PSEE7]. 
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The purpose of the present work is to extend the theory of regularity structures to 
multicomponent systems of the form 


d t u = A x u + F(u, v) + £ , 

dtv = G(u,v). (1.3) 

The main result is the proof of existence, in space dimensions d = 2 and 3, of local solutions 
to the system (11,31) when F is a cubic polynomial and G is linear. This includes the case 
of the standard FitzHugh-Nagumo model, but also other equations such as the Koper 
model |2Bj with diffusion in the fast component, which features vectorial variables v. We 
will mostly focus on the case d = 3, first because this is the physically relevant case in 
many applications, but also because it is the technically more challenging case (the Allen- 
Cahn equation is known not to be renormalisable for d ^ 4, so the same will hold a fortiori 
for multicomponent systems with cubic nonlinearities). 

The main technical difficulty that has to be overcome to obtain these results comes 
from the fact that the semigroup associated with the second equation in (11.31) is not at 
all regularising in space. Though one expects that this loss of regularisation is somehow 
compensated by the fact that no singular noise term acts on the equation for v, the 
general theory in m cannot be applied directly, because it uses in an essential way the 
assumption that the heat kernel is smooth everywhere except at the origin. Therefore, 
we have to extend the regularity structure constructed for the Allen-Cahn equation with 
new abstract symbols representing integration with respect to a singular kernel, and to 
prove suitable bounds involving these symbols. Furthermore, one has to analyse which 
role new symbols play in the renormalisation. We find that certain terms involving the 
singular kernel do not have to be renormalised, while on the other hand, general cubic 
nonlinearities yield renormalisation terms which do not appear in the Allen-Cahn case. 

The remainder of this work is organised as follows. Section [2] contains all main local 
existence results. Section [3] gives a summary of the most important aspects of the theory of 
regularity structures contained in m, illustrated in the case of the Allen-Cahn equation. 
In Section lU we present our results allowing to represent the operation of integration with 
respect to a singular kernel. Section [5] contains the fixed-point argument proving local 
existence and uniqueness of solutions for the SPDE with mollified noise, and Section [6] 
deals with the renormalisation procedure which is necessary to pass to the limit of vanishing 
mollification. Finally, in Section [7] we complete the proofs of the main results. 

Notations: We write |.t| to denote either the absolute value of x £ M or the I 11 -norm of 
x £ M. d , while ||x|| denotes the Euclidean norm of x £ M. d . If a, b £ R we use a A b := 
min{a, b} and a V b := max{a, b}. If /, g are two real-valued functions depending on small 
parameters e, 6,. .. (which will be clear from the context), we write / < g to indicate that 
there exists a constant C > 0 such that / ^ Cg holds uniformly in the small parameters. 
We use the notation f x g to indicate that we have both / < g and g < /. 

Acknowledgements: N.B. would like to thank the Institute for Analysis and Scientific 
Computing at the Technical University in Vienna for kind hospitality and financial sup¬ 
port. C.K. would like to thank the Austrian Academy of Sciences (OAW) for support 
via an APART fellowship. C.K. also acknowledges support of the European Commission 
(EC/REA) via a Marie-Curie International Re-integration Grant. 
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2 Results 


2.1 The standard FitzHugh—Nagumo equation 

We start by considering the particular case of the standard FitzHugh-Nagumo SPDE, 
given by 


dtu = A x u + u — u 3 + v + f , 

dtv = a\u + a2V. ( 2 - 1 ) 

Here A x = J2j=i ®xj is the usual Laplacian, u and v are functions of time t ^ 0 and space 
x £ T d (the torus in dimension d = 2 or d = 3), £ stands for space-time white noise on 
R x T d , and 01,02 a real parameters. As such, the system ( 12 . 11 ) is not well-posed, and a 
renormalisation procedure is required to define a notion of solution. To do this, we first 
choose a rescaled modifier 

&(*.*) = ^2 e (j?'f) ’ (2 ' 2) 

where g : R rf+1 —>• M is a smooth compactly supported function of integral 1. We set 
= Qe * where the star stands for space-time convolution, and consider the sequence 
of equations 


d t u £ = A x u £ + [1 + C(e)] u £ - ( u E f + v £ + f £ , 

dtv £ = a\u £ + a, 2 V £ , (2.3) 

where C{e) G R. Then our first result is the following, which is close in spirit to [19, 
Thm. 1.15]. 

Theorem 2.1 (Standard FitzHugh-Nagumo SPDE). Assume uq belongs to the Holder 
space C 71 for some rj > — |, and uo belongs to the Holder space C 7 for some 7 > 1. Then 
there exists a choice of the constant C(e) such that the system (12.31) with initial condition 
(uq,vo) admits a sequence of local solutions ( u £ ,v £ ), converging in probability to a limit 
(u,v) as e —> 0. This limit is independent of the choice of mollifier g. 

The proof of this theorem and the two next ones is given in Section [7| For the precise 
definition of Holder spaces C v with negative index, we refer to Definition 13.31 below. By 
local solution we mean that for any cut-off L > 0, the solution is defined up to the random 
time when the C^-norm of ( u £ ,v £ ) first reaches L. In other words, we cannot exclude the 
possibility of finite-time blow-up. 

The constant C(s) admits an explicit expression in terms of the heat kerneQ 

G " i,(M) = li^ e " w2/4 ' 1 i M '' (2 ' 4) 

In dimension d = 3, C(e) is the same as for the dynamical <f >3 model considered in [19]. 
Namely, setting = G^ * g e , one has 

(7(e) = 3 f G^(z) 2 dz- 18 f G^{z)Q £ 0 (z) 2 dz , (2.5) 

Jr 4 c Jr 4 

1 Below, the renormalisation constants will rather be defined in terms of a truncated version of the heat 

kernel. This does however not affect their singular parts (terms of order e _1 and log(e -1 )). 
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/o\ /o\ 

where Qq(z) = f Ge (zi)Gl (zi — z) dz\. In particular, it has been known since the works 
of Feldman and Osterwalder mm that 

C(e) = ^ + C 2 log(e” 1 ) + C 3 ( 2 . 6 ) 

for some constants C\,C 2 , C 3 € M. The constants C) and C 3 depend on the choice of 
modifier g, while C 2 is independent of g (cf. [T51 Rem. 6.2]). 

In dimension d = 2, one can show (cf. [23l Rem. 2.14]) that the renormalisation 
constant is given by 

COO = 3/ G?\z) 2 dz = ^-log(e- 1 ) + C 3 , (2.7) 

J R 3 47r 

for some constant C 3 G M, depending again on the choice of g. 

The proof of Theorem 12.11 will also provide some information on the structure of the 
solutions. Indeed, we always have 

u £ {t, x) = Xe(t, x) + ip e {t, x) , (2.8) 

where Xe = G ^ *£ £ is the stochastic convolution of heat kernel and noise, and ip e converges 
to a function (as opposed to a distribution). In other words, the only singular term in the 
limit e —> 0 is given by lim e _>.o \e = * £, which is independent of the nonlinear term 

in the equation. The function v £ (t,x) has the same structure, as it can be represented in 
terms of u e by solving a linear inhomogeneous equation. 

2.2 More general cubic nonlinearities 

We can now extend the above results to a more general class of systems, of the form 

d t u = A x u + F(u, v) + £ , 

dtv = a\u + 02 U , (2.9) 

where F is a cubic polynomial of the form 

F(u, v ) = a\u + a^v + /3±u 2 + fauv + fov 2 + 71 n 3 + X 2 U 2 v + 73 uv 2 + 74 U 3 . (2.10) 

In that case, the renormalised equations take the form 

d t u £ = A x u £ + [F(u £ ,v £ ) + Co(e) + Ci(s)u £ + C 2 (e)u £ ] + , 

dtv s = a\u £ + a 2 v e , ( 2 . 11 ) 

where Co(e),Ci(e),C 2 (e) € K, and we have the following result: 

Theorem 2.2 (General cubic nonlinearities). Assume uq and vq satisfy the same as¬ 
sumptions as in Theorem \2.1[ Assume further that either d = 2, or d = 3 and 72 = 0. 
Then there exists a choice of constants C'o(e), C'i(e) and C 2 (e) such that the system (12.111) 
with initial condition (uq,vo) admits a sequence of local solutions (u £ ,v £ ), converging in 
probability to a limit (it, v) as e 0. This limit is independent of the choice of mollifier g. 
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Note that in dimension d = 3, we assume that F contains no term of the form r y 2 u 2 v . 
This is because the method we use does not yield a simple form for the renormalised 
equation if such a term is present (see Section El, which is an indirect consequence of 
the fact that the equation for v is not regularising in space (cf. Section 14.211 . It is not clear 
at this stage whether this is just a technical artefact, or whether it has a deeper meaning. 

The renormalisation constants can again be computed explicitly to leading order. They 
are given by 

C 0 (e) = -±p 1 C(e), 

Ci(e) = - 7 iC(e), (2.12) 

C 2 (e) = -p 2 C(e), 

where C(e) is given by 

C(e) = 3 [ Gf\zf &z + 1%^ f G^Xz)Ql(z) 2 &z (2.13) 

J R 4 J R 4 

for d = 3 and by (12.71) for d = 2. Note in particular that these constants depend only on 
the coefficients of u 2 , u 2 v and u 3 in F. 


2.3 Vectorial gating variables 


Another generalisation of interest is to systems with vectorial variables v, of the form 

d t u = A x u + F(u, v) + £ , 

d t v = uAi + A 2 v , ( 2 -14) 


where v(t,x) takes values in M n . Here A\ G is a constant vector, and A 2 G K nxn is a 
square matrix, while F is again a cubic polynomial in u and the u*. This allows for instance 
to consider the Koper model [26] with spatial diffusion and space-time white noise in the 
fast variable; for the stochastic Koper model without the Laplacian representing spatial 
diffusion see [5]. The reaction terms of one version of the Koper model can be written in 
the form (12.141) with 

F(u, v) = 3u + vi - u 3 , A, = (j^') i ^2 = ^ ^ , (2.15) 

where k,e± Gl are model parameters. 

The natural candidate for the renormalised system associated to (12.141) is given by 


dtu e = A x u £ + 


F{u e , v £ ) + Co(e) + Ci{e)u £ + ^ ] C 2 ^(e)i 


i=l 


+r, 


d t v £ = u £ Ai + A 2 v £ , 


(2.16) 


for constants Cq (e), C'i(e), C 2 ^ (e) G M. Indeed, we have the following result. 

Theorem 2.3 (Vectorial variables v). Assume uq and the components of vq satisfy the 
same assumptions as uq and vq in Theorem \2.H Assume further that either d = 2, or 
d = 3 and F(u,v) has no terms in u 2 Vi. Then there exists a choice of constants Cb(e), 
Ci(e), C 2t i(e) such that the system ( 12 . 161 ) with initial condition (uo,vq) admits a sequence 
of local solutions (u £ ,v £ ), converging in probability to a limit (u,v) as £ —> 0. This limit 
is independent of the choice of modifier g. 
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The renormalisation constants can again be expressed in terms of the coefficients of 
the initial equation (|2.14p . Writing 

n 

F(u , v) = aiu + Piu 2 + 71 n 3 + ^ [a 2 ,i ?7 + & 2 ,iUVi + 72 ,iU 2 Vi] + R(u, v ) (2.17) 

i= 1 

with u)| ^ C||u|| 2 , we have 

Co(e) = -\(hC{e) , 

Ci(e) = -71 C{e) , (2.18) 

C 2 ,i{e) = -- 72 ,iC'(e) , z = 1 ,..., n , 

where C(e) is again the constant defined in (J2.13I) for d = 3 and in dSZD for d = 2 . 
Note that in the particular case of the Koper model, C\{e) is equal to the constant C(e) 
obtained for the FitzHugh-Nagumo equation, while all other constants vanish. This is not 
surprising, since both models have the same nonlinearity. 

Remark 2.4. The notion of local subcriticality given in m Assumption 8.3] suggests that 
the class of renormalisable SPDE-ODE models of the form (jl.3l) is larger than the one 
considered here. More precisely, in dimension d = 3, one would expect models with quartic 
F and linear G to be renormalisable, as well as models with F quartic in u and quadratic 
in v and G quadratic in u and linear in v. The reason we did not include them in our 
analysis is a technical one: quartic F would produce solutions whose singular component 
is not homogeneous in space and time, as in (12.81) . The way in which we lift the singular 
kernel to the regularity structure is not able to deal with such situations (see Section 14.21 
below). We plan to further investigate this issue in future work. 0 

3 Regularity structure for the Allen—Cahn equation 

This section serves the double purpose of giving a very brief account of the theory of 
regularity structures contained in pj2J, and of describing a regularity structure for the 
Allen-Cahn equation 

dtu = A x u + u — u 3 + £ , (3-1) 

where A x = J2 C j=i^xj is the usual Laplacian, £ = £(t,x) denotes space-time white noise 
on R x T d where denotes the d-dimensional torus with d = 2 or d = 3, and we seek a 
solution u : [0, T] x —> M for a given initial condition uq = tt(0, x ) in a suitable function 
space. The regularity structure for (I3.ip will serve as our starting point to build a larger 
structure allowing to represent the FitzHugh-Nagumo equation. 

By Duhamel’s principle, one possible solution concept for (13.11) is to consider the fol¬ 
lowing integral equation 

u t = [ S(t-s)[u s -ul + £ s ]ds + S(t)u 0 , (3.2) 

Jo 

where S(t) = e tAa: denotes the semigroup of the heat equation compatible with the bound¬ 
ary conditions and u s ■= u(s, ■), := £(s, •). The purpose of a regularity structure is to 

provide an abstract space in which one can construct a fixed point of (13.21) when £ is 
replaced by a mollified version £ £ . Then the idea is to take the limit e —> 0, and to project 
the solution to a distribution on [0, T] x T rf , where “distribution” is understood in the 
sense of a “generalized function” representing a sample path. 
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3.1 Regularity structures 

Definition 3.1 ( [T9l Def. 2.1]). A regularity structure is a triple ( A , T, Q) consisting of 
(Rl) an index set A C M, containing 0, which is bounded from below and locally finite; 

(R2) a model space T, which is a graded vector space T = 0 aeA where each T a is a 
Banach space; the space Tq is isomorphic to R and its unit is denoted 1; 

(R3) a structure group Q of linear operators acting on T, such that 

Tt ~ t € © T /3 = : T a (3-3) 

P <a 

holds for every T £ Q, every a £ A and every r £ T a ; furthermore, T1 = 1 for 
every T € Q. 

Example 3.2. A simple but important example of regularity structure is the polynomial 
regularity structure for d variables. In that case, A = No is the set of non-negative integers. 
For £ £ No, Tg is the space of homogeneous polynomials in d variables of degree t. It is 
spanned by the monomials X k = X kl ... X kd for which |fc| = ki + • • • + kd = £■ Finally, 
the structure group Q is defined by 

F h (X k ) = {X - h) k , h £ • (3.4) 

This group is isomorphic to R d , and one easily sees that satisfies the requirement (I3.3p . 
The interpretation of Q is that it allows to convert a Taylor expansion around a point 
x £ R rf into the expansion around another point x + h. ♦ 

We will henceforth denote by T the polynomial regularity structure with d+ 1 variables 
Xq, ..., Xd, where Xq represents the time variable. Since the linear part of the Allen-Cahn 
equation is given by a parabolic operator, it turns out to be useful to make the time variable 
“count double”. This is done by defining the parabolic scaling 

S = (2,l,...,l)€N d+1 , (3.5) 

acting on M rf+1 and declaring that a monomial X k has homogeneity |A' fc | s = |/c| s where 
the scaled degree is defined by \k\ 5 = 2k$ + ^T =1 ki. 

The regularity structure of the Allen-Cahn equation (|3.3I) in M. d is built by enlarging 
T, i.e., by adding new symbols other than the polynomial symbols X k to T (see [23]). The 
noise is represented by a symbol 5. In order to account for the fact that space-time white 
noise has Holder regularity a for any a < — (d + 2)/2, we set 

d + 2 

a 0 = --- k (3.6) 

where k > 0 will be chosen sufficiently small in the sequel, and declare that H has homo¬ 
geneity |H| S = a 0 . 

The set of symbols is equipped by a product, which by definition is commutative 
and associative with unit 1. The product of two elements r, a £ T has homogeneity 
\to\ s := |r| s + |cr| s . Furthermore, integration against the heat kernel is represented by a 
map X : T —>• T, which by definition satisfies |X(r)| s := |t| s + 2, in order to account for 
the regularizing effect of the heat kernel. 

Let T be the set containing all possible products of symbols 1,5, Xj and their images 
by X, and denote by XL = span(X”) the vector space spanned by all these symbols. This is 



T 

Symb 

|t|. 

d = 3 

d = 2 

A(t) 

77 

77 

OiQ 


—2 — K 

H <g> 1 

Z(H ) 3 

•v 

3ckq + 6 

-3« 

0 — 3k 


m 2 

V 

2 «o + 4 

-1 - 2 k 

0 — 2 k 

V <g> 1 

Z(Z(Z) 3 )Z(Z) 2 

•Y- 

5ao + 12 

-i -5 k 

2 - 5k 

■Y* 8>1 + V ®77(^) 

1(Z) 

T 

ot 0 + 2 

“3 ~ K 

0 — K 

T <g> 1 

Z(Z{Zf)Z{Z) 

Y- 

4ao + 10 

0-4 k 

2 — 4k 

Y* <g> 1 + T ®J(W) 

Z{Z{Zf)Z(Zf 

■Y- 

4ao + 10 

0 — 4k 

2-4k 

•Y* 8 > 1 + V ® 77 (v) 

X(H) 2 W 

VI, 

2 ckq + 5 

0 — 2 k 

1 - 2 k 

V ® Xi + Vlj ® 1 

1 

1 

0 

0 

0 

1 ( 8 ) 1 

l(l(Zf) 

T 

3ao + 8 

3 ~ 3 k 

2 — 3k 

Y < 8)1 + 1 ® J(^) 

Z(Z(Zf)Z(Z) 

Y* 

3ao + 8 

3 — 3k 

2 — 3k 

Y* ® 1 + T <g>J(V) 

Z(Z{Z))Z(Zf 

•X* 

3ao + 8 

1 -3k 

2 — 3k 

*Y* 8) 1 + 8> 77( T )+ 

VXi®Ji( T) + V ®X^(T) 

Z(Z(H) 2 ) 

Y 

2oq + 6 

1 - 2k 

2 - 2k 

Y ®1 + 1 8>J(V) 

Z{Z{Z))Z(Z) 

U 

2ao + 6 

1 - 2k 

2-2 k 

L* ®1+ T 8> t 7(T)+ 

iXi®ji(\) + i®Xiji(\) 

Xi 

Xi 

1 

1 

1 

Xi®l + l®Xi 

Z(Z(Z)) 


op T 4 

| -K 

2 — K 

® 1 + 1 ® i7(T)+ 

Xi®Ji{\) + l®XiJi(]) 


Table 1. Elements of Tp of lowest homogeneity for the Allen-Cahn equation. They have 
been ordered by increasing homogeneity for the case d = 3. The expressions for A(r) are 
shown for the case d = 3. We have written J instead of Jq and Ji instead of J ei , where 
the d are canonical basis vectors of Z+ +1 . Summation over the index i is understood. 


an infinite-dimensional space but in practice only a finite-dimensional subspace of T-L will 
be needed. Let U be the smallest set containing 1, Xi and 21(H), and such that 

Ti,T 2 ,T 3 eU =7 1(tit 2 t 3 )£U. (3.7) 

We then set 

7> = {H} U {tit 2 t 3 : Ti£U} , (3.8) 

and define the model space as being the vector space T = span(J 7 p) spanned by Zp. The 
index set is then defined as 

A = {|r| s : t € Zpj • (3.9) 

The model space T admits a natural grading 

T = ®T 7 (3.10) 

7 eA 

obtained by letting T 7 be the vector space spanned by elements of homogeneity 7 . We 
equip each T 7 with a norm ||-|| 7 ; the choice of norm is irrelevant since pfQi Lemma 8.10] 
shows the (nontrivial) fact that all T 7 are finite-dimensional. 
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Table [T] shows the elements in Tp of lowest homogeneity, using a graphical represen¬ 
tation introduced by Hairer. Each symbol S is denoted by a dot, and integration with 
respect to I is denoted by a vertical line pointing downwards. For instance, X(E) = T. 
Multiplication of symbols is represented by joining them at the base, so for instance 
1(E) 2 =V. 

Finally, we have to extend the definition of the structure group Q, which is the most 
involved part of the construction. The first step is to introduce a set J-+ C F, which 
contains all formal expressions of the type 

xk U Jk^j, (3.11) 

3 

where tj G F and the multiindices kj are such that |tj[ s + 2 — \kj\ s > 0. By definition, an 
expression of the form (13.1111 has homogeneity \k\ s + J2j(\ T j\s + 2 — \kj\ s ), which is always 
strictly positive, except for the element 1, which has homogeneity 0. The set Fp C F+ is 
defined similarly, but with Tj G Fp. We set T-L+ = span(J r + ) and T + = span(J 7 i )i). 

The last ingredient is given by a bilinear map A : T~i —> H <8 W.+ defined by setting 

Al = l(g)l, AX* = Xi® 1 + 1®X* , AS = S <8 1 , (3.12) 

and extended inductively to all of 'H by the rules 

A(t<t) = (At)(Acj) , 

_ j^£ 

A {It) = (X <8 Id)A(r) + Y] —jr <8 — rJi+ m r , (3.13) 

Z —' t\ ml 

£,m 

where Id denotes the identity map and we impose that 

J k T = 0 if |fe| s > \t\ s + 2 , (3.14) 

while J^t is a new formal symbol otherwise. Table Q] also shows A(r) for the first few 
elements of Fp. 

Denote by Tip the dual of H+, that is, the set of linear maps h : T-L+ M, that 

we will write as r G (h, r). Let Q+ be the set of grouplike elements g G %*, that is, 
those satisfying (g,TC r) = (g,T)(g,cr) for all r, cr G Then Q is obtained by identifying 
elements of Q± acting the same way on T + . The action of Q on T is defined by 

(g,r) i-G V g T = (Id( 8 ff)Ar . (3.15) 

The fact that (A,T,Q) constructed in this way is indeed a regularity structure is proved 
in [191 Th ru . 8.24], 

3.2 Models and reconstruction theorem 

A regularity structure can be used to represent an SPDE as a fixed-point problem in an 
abstract space. One also has to provide a connection between the abstract space and 
the “physical” space the solution lives in. First, we need to introduce a version of the 
classical Holder spaces C a on In order to take the parabolic scaling into account, 

the Euclidean norm is replaced by 

\\(t,x)\\ s = |t | 1/2 + 5^1*4 - (3.16) 

i =1 
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Denote by (1, •) the element of the dual of T defined by (1,1) = 1 and (l,r) = 0 for all 
r G ® 7 _£ 0 T 7 . If a > 0, the space C“ is defined as the set of all functions ip : R rf+1 —> R 
such that there exists a function <p : R rf+1 —>■ T~ such that (1 ,(p{z)) = ip(z) for every 
z = (t. x ) G R rf+1 and such that for every compact set C R rf+1 the estimate 

\<p(z + h)-r h <p(z)\p <\h\“~P (3.17) 

holds uniformly over (5 < a, \h\ s ^ 1 and z G Jl (cf [19l Def. 2.14]). Lemma 2.12 in [19] 
shows that in the case of the Euclidean scaling s = (1,..., 1), this definition coincides with 
the usual definition of Holder spaces. 

If a < 0, Cg can be defined in a natural way as a subset of Schwarz distributions 
(S^R'L -1 ) as follows. Given r G N, denote by B r s 0 the space of functions of class C r which 
are supported in the set {z G R rf+1 : ||z|| s ^ 1}. 

Definition 3.3 (j 49l Def. 3.7]). Assume a < 0 and let r = — |_aj. A Schwarz distribution 
v G <S'(R d+1 ) belongs to if it belongs to the dual of Cq and for every compact set 
A. C R rf_,_1 ; there exists a constant C such that 

(v.Si.rt) ^ C5° (3.18) 

holds for all rj G B r s 0 with \\r]\\cr ^ 1, all 5 G (0,1] and all z G Jt. Here 

: = 5~ [d+2) ri(5~ 2 {i-t),5- l {x-x)) . (3.19) 

Next we introduce the notion of a model associated with a regularity structure (H, T, Q). 

Definition 3.4 ( fl9l Def. 2.17]). A model for a regularity structure (A , T, Q ) with scaling s 
is a pair (n,r), defined by a collection {n^ : T —> 5 / (R rf+ 1 )} zgR d+i of continuous linear 
maps and a map T : R rf+1 x R^ -1-1 —>• Q with the following properties. 

(Ml) T zz = Id is the identity of Q and T zz >T z i z n = F zz » for all z,z',z" G R d+1 . 

(M2) U z , = H Z F ZZ , for all z , z' G R d+1 . 

(MS) For any 7 G R and any compact set C R d+1 , one has 
||n|| 7; * := sup sup sup sup sup 

ze& a<7 reT a 0<5^1 

(Mf) For any 7 G R and any compact set & C R rf+1 , one has 

ii-pii \\ r zz'T\\l3 / 001 \ 

||1 || 7; ^ := sup sup sup sup - a _3 < 00 • (3.21) 

z,z'£R a <7 f3<a reT a ||r ||\\z — Z r \\s 

In the case of the Allen-Cahn equation with mollified noise := £ * g £ , where g(t,x ) 
is a mollifier and g £ (t,x ) = e~^ d+2 ^ g(t / e 2 , x/e), a canonical way of building a model 
Z £ = (Hy F £ ) proceeds as follows. First, we can just define 

(H £ E)(f) = e(z) , 

(H £ z X k )(z) = (z-z) k , 

(H £ z tt)(z) = (I W z t)(z){U £ z t)(z ) Vr, t G T . (3.22) 


\(IL z T,Sl z T})\ 


< 00 . 


(3.20) 
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To also include the integration map in the fixed point problem let G = G(t, x ) be the heat 
kernel defined as the fundamental solution associated to dtu = A x u. We split the heat 
kernel G as 

G = R + K , (3.23) 

where K : (R d+1 \{0}) —> R is a singular part satisfying specific algebraic properties, while 
R : R d+1 -> R is a smooth part. The properties of K are: 

• K is supported in {|x| 2 + \t\ ^ 1} where |x| = Y^j=i\ x j\'i 

• K(t,x) = 0 for t ^ 0 and K(t, —x) = K(t,x ) for all (t, x); 

• we have 

K &= ]^\dj 2 e-"" l|2/4 * for |x| 2 + |tK \ (3-24) 

and K(t,x ) is smooth for |x| 2 + \t\ > 

• furthermore, one has vanishing moments 


/ 

J R d +! 


K(t, x)P(t, x) dx d t = 0 


(3.25) 


for all polynomials P of parabolic degree less or equal some fixed £ ^ 2. 

Lemma 5.5 in m shows that such a splitting indeed exists and also the vanishing 
moments condition holds for K ; both parts K,R then satisfy a number of derivative 
bounds. Furthermore, it follows from Lemma 5.5 in |19| that there exists a decomposition 


K(z - z') = y, K n (z - z') 
n^O 


(3.26) 


where n G No and suitable derivative bounds and vanishing moment conditions hold for 
each of the kernels K n . With this construction, the abstract integration map is represented 
by the Taylor-series-like expression 

ffil r)(z) = f K(z - z')(n% T )(&z') + Y, { -^ —UlJir) , (3.27) 

J i 

where the linear forms /5 G Q C T~L\ are constructed as follows: 


(f z ,Xi) = -Zi , 

(. fz, Tf ) = if!’ T ) Vr.fGT, 

(/I, JlT) = ~j D e K(z - z)(m T )(dz) . (3.28) 

Note that in the more general kernel case K = I\(z, z) one has to replace D^ by D\ and 
the subscript indicates derivative with respect to the first argument. The second term 
in the definition (13.2711 of II e z Xr may seem somewhat strange, and will be motivated in 
the next section. Note that it should really be interpreted by applying both sides to test 
functions, i.e., for all smooth compactly supported functions ip and all r G T a one requires 

(n^T, </>} = £/ V^'xnir, K% ,) dz 1 (3.29) 
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where K%. zz , is defined via 

K. zz ,(z") = K n {z'- z") - Y, {Z '~ Z)k D k K n {z-z"). (3.30) 

\ k \ s < a-\-2 

The group elements T e z % are then simply defined by 

r % = (FI)- 1 FI where F £ z = T /f (3.31) 

(recall that FI is invertible because f £ z is grouplike). In this way, the algebraic property 
(Ml) of Definition 13.41 is automatically satisfied. Regarding the algebraic property (M2), 
we can note the following relations. Using Sweedler’s notation [35], we write At = t*' 1 ! <S> 
r 2 \ although At is usually a sum of such terms. Let further j zz denote the element of 
the structure group such that T £ zz = pc . Then we have by (13.151) 

mrL-r = II^(Id 07 ^)At = n. (3.32) 

Property (M2) thus amounts to the relation 

n i r = n ^ (1) ( r ^ r(2) )> ( 3 - 33 ) 

which provides some intuition for the meaning of A. It is easy to check that this relation 
holds for elements r of the polynomial regularity structure. In the general case, the fact 
that Z £ = (IF,r £ ) is indeed a model is a nontrivial fact, proved in |I9| , Prop. 8.27]. 

Remark 3.5. It is very important to realise that the canonical model just built is not the 
only possible model for a given regularity structure. This freedom in the choice of model 
will be used when introducing the renormalisation procedure. All models, however, will 
share many properties with the canonical model. The only rule that will be modified is 
the product rule II j(tt) = n £ (r)II|(r). 

We can now introduce the spaces D 7 , which play an analogous role as the C“ on the 
level of the regularity structure. 

Definition 3.6 ( |19 . Def. 3.1]). Let 7 £ R. Given a model Z = (II, T), the space X > 7 = 
T>' y (Z) consists of all functions f : M rf+1 — > T~ such that for every compact set R C M rf+1 
one has 

II/IIIt;^ : = sup sup||/(^)|| / a + sup sup < 00 • ( 3 - 34 ) 

zGJt /3<7 z,z£fi /3<7 \\Z — Z\\l 

Wz-zWs^l 

In the particular case of the polynomial regularity structure, it is again quite straight¬ 
forward to check that the requirement (13.341) is equivalent to / being the Taylor expansion 
of an element of the Holder space C]. The spaces D 7 depend on the model via T, but not 
on n, cf |19( Remark 3.4]. In order to compare elements of X> 7 for different models, it is 
useful to introduce 


MHI/-/IU + 


sup 

z,z£& 


sup ■ 

/3<7 


\\f(z)-f(z)-r gg f(z) + r eS f(z) 


z — 


~ll 7-/1 


(3.35) 
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(which is in general not a function of / — /), where 


11/ - /ll 7 ;.« = SU P SU P \\f( z ) ~ f( z )b ■ ( 3 - 36 ) 

z&R. / 3<7 

The central result allowing to link elements in T> 7 and in C“ is the reconstruction 
theorem US Thm 3.10]. It states that if a* = min .A and given r > |a*| there exists, for 
any 7 £ M, a continuous linear map TZ : T> 7 —>■ C"* such that 

\(7Zf-n z f(z),S s s z r,)\ ^C<5 7 ||n|| 7;I ||/| 7;I (3.37) 

holds uniformly over all test functions r) £ B s r 0 , all 6 £ (0, 1 ], all / £ V 7 and all z £ M. 
The constant C > 0 depends only on 7 and on the regularity structure, and the set ll is 
the 1-fattening of M (i.e., the points at distance at most 1 from J4). If 7 > 0, then (13.371) 
defines TZ uniquely. The reconstruction theorem also provides bounds on the dependence 
of TZ on the model (II, T). Heuristically, (13.371) states that 7 Zf locally looks like U z f(z) 
near any point 2 £ K d+1 , up to terms of order 7 . 


3.3 Lifting the convolution maps 

We can now lift the operation of convolution with the heat kernel to the space T> 7 (Z). 
The requirement (13.251) is for compatibility with the condition that X should define an 
abstract integration map of order 2 in the sense of |19| Def. 5.7], namely 

• X : T a —> T a .|_2 for every a £ A] 

• Xt = 0 for every r in T; 

• XTt — TXt £ T for every r £ T and every T £ Q. 

The central result is given by the so-called multilevel Schauder estimates [ 191 Thm 5.12], 
which state in particular that for all 7 £ M such that 7 + 2 0 N, there exists a map 
/C 7 : T > 7 —> £> 7+2 such that 

IZJC^f = K * XZf ( 3 . 38 ) 

holds for all / £ X) 1 . In other words, the following diagram commutes: 


D 7 


/C, 


£>7+2 


TZ 


TZ 


cr 


K* 




The map /C 7 has the following expression. For any / £ D 7 , 


(/C 7 /)(z) = Xf{z) + J(z)f(z) + (A/- 7 /)(s) , 


where for each r £ T a , 

J(z)t = 

W){z) = 


E 

x k 

Id 

[ D k K(z 

J K d +! 



E 

X k 

IT, 

[ D k K(z 
J R d +! 

|fc| s <7+2 




14 


z)(JI z t){&z) , 
z)(TZf -U z f(z))(dz) . 


(3.39) 

(3.40) 

(3.41) 











Note that the last two operators have values in T, the polynomial part of the regularity 
structure, and that the only nonlocal operator is A7y. 

The role of the operators J and A /" 7 is to ensure that /C 7 / has the properties required 
to belong to T> 7+2 , which would not be the case if one simply sets /C 7 / = Tf. The maps 
J are related to the coefficients Jk appearing in (13.111) : in fact, the Jk T j play the role 
of placeholders for the J[z)t. 

As for the smooth part R of the heat kernel, it can be lifted as in [ 191 (7.7)]. Namely, 
with a smooth kernel R we associate maps i ? 7 : —> V 1 given by 

(Ryv)(z)= Y Ti" / D k R(z-z)v(z)dz:= Y D kR ( z - 0) ■ ( 3 - 42 ) 

|fc| s <7 ' ® + |fc|s<7 

It follows from m Prop 3.28] that for v € one has 

{TZR 1 v)(z) = (v,R(z - •)) , 


and thus 

IZR^f = R*Kf . 

In other words, the following diagram commutes: 


RGR 

W - - -* TP 



(3.43) 


Assume that for an appropriate choice of 7 , one can find an element U £ D 7 (Z e ) 
satisfying the fixed-point equation 

U = (/C 7 + R^R+tZ + U — U 3 ) + Gu 0 , (3.44) 

for some 7^7 — 2 , where R + (t,x) = lp>o} an4 Guq denotes a suitable lift of the 
convolution in space of heat kernel and initial condition. Applying the reconstruction 
operator 1Z to both sides of this equation, one can show that u = 1ZU satisfies 

u = (K + R) * (. R + [u - u 3 + £ £ ]) + Gu 0 , (3.45) 

which is equivalent to (13.21) ; in (13.45[) we also use Guq as a notation for the usual convolu¬ 
tion in space of heat kernel and initial condition. This is basically the strategy implemented 
in Theorem 7.8 and Section 9.4 of |19| . except that one has to deal with two additional 
technical difficulties. The first one is that due to the singular behaviour of the heat kernel 
as time goes to 0, the definition of the space D 1 has to be modified. We will apply this 
modification to our case in Section I4~3l The second difficulty is that although a fixed point 
exists for every mollification parameter e > 0, it cannot converge as e —> 0. This is why a 
renormalisation procedure is needed, which we will adapt to our case in Section [ 6 j 
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4 Extension of the regularity structure 


Our aim is now to extend the regularity structure built for the Allen-Cahn equation, in 
order to allow to represent a family of coupled SPDEs-ODEs of the form 


d t u = A x u + F(u, v) + , 

d t v = uA\ + A 2 v , (4.1) 


where F(u, v ) is a cubic polynomial, and A± £ R ra , A 2 £ M nxn are either scalars as for the 
classical FitzHugh-Nagumo case with n = 1, or a vector and a matrix if v has multiple 
components. Duhamel’s formula allows us to represent (mild) solutions of (|4.1I) as 


u t = [ S(t- s)[£f + F(u s ,v s )] ds + S(t)uo 

Jo 

v t = / u s Q(t - s) ds + e t4z vq , 

Jo 


(4.2) 


where Q(t) := e tA2 A\. We thus have to lift to the regularity structure the operation of 
time-integration with respect to Q , which has no smoothing effect in space. 

In the case where v has values in R n with n > 1, the kernel Q(t) is in fact a vector 
of dimension n. In what follows, we will mainly deal with the case of scalar Q , since the 
generalisation to the vectorial case is rather straightforward. 


4.1 Extension theorem 

Let us fix a finite time horizon T. Then we can always assume that the kernel Q satisfies 
the following properties: 

• Q(t) is supported on [0,2T] and smooth for t > 0; 

• Q(t) = e tA2 A\ for all t £ [0,T], 

The reason why this is allowed is that we will be interested in showing existence of solutions 
on a sufficiently small interval [0, T], so that the behaviour of Q outside this interval is 
not going to matter. Also, since Q is bounded, it will not be necessary to decompose it 
as a sum of Q n concentrated in sets of radius 2~ n , as in the case of the heat kernel in [191 
Section 5], cf (13.2611 . 

In order to be able to represent the time-integration map, we will have to extend our 
regularity structure and the associated model. We do this by adding to T new symbols 
denoted £(r), r £ T \ T, which we represent by an open blue dot. Thus for instance we 
write £(Z(S)) = £(T) = ?. In practice, we will only need to apply £ to elements r of 
homogeneity |r| s £ (—2,0). We thus set V = {r £ T: — 2 < |r| 5 < 0} (V is called a 
sector of T). If r 0 V, we simply set £{r) = 0. We also postulate that for r £ V, £{t) 
has the same homogeneity as r, and we choose the norm on the vector space generated 
by the new symbols in such a way that 

||£(t)|| q = ||t|U Va £ E . (4.3) 

The operator £ defined in this way is an abstract integration map of order 0 on V in the 
sense of [19l Def. 5.7], i.e., 

• £ : V fl T a —>• T a for every a € A; 

• £t = 0 for every r £ V fl T; 
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• £Tt — T£t £ T for every r £ V and every T £ Q. 

The first two properties are obvious by definition of £t, in particular we have V (IT = 0. 
Regarding the third property, we first extend the structure group by setting 

A(£t) = {£ <g> Id)Ar . (4.4) 

The third required property of £ is indeed satisfied since 

T g £r = (Id <g>g) A(£r) = (£ <g> g) At = £T g r , (4.5) 

so that in fact T^t — £T g r = 0 for all r £ T and all g £ Q. 

Since £ is not regularity-increasing, there can be in principle infinitely many symbols of 
given homogeneity. The trick, however, will be to build the fixed-point map in such a way 
that only finitely many new symbols are needed. In practice, it will turn out that the only 
required new symbols are ?, and those obtained by applying T to existing symbols and/or 
multiplying them. However, since it is of independent interest, we are going to describe 
the extension procedure in a more abstract, inductive way. Given a subset W C V, on 
which the canonical model Z £ = (n e ,T e ) is defined, we want to extend the model to a 
larger set W = W U {£{j) : r £ W}. We can then apply the usual extension theorem [19, 
Thm 5.14] to extend the model to W U {X(r): r £ W}, and so on as often as needed. 

The inductive step from W to W goes as follows. Assume that Z e already satisfies 
Definition E3] on a regularity structure (A, W, Q). We define the extended model by setting 

( n t,aA)(f>®) = [_ Q(i-s)(U £ tx T)(s,x)ds (4.6) 

Jt-2T 

for all r £ W. Writing as before F £ = Tye, the new group elements are defined by setting 

T £ zz ,(£t) = (. F £ z )~ l F £ A£r) (4.7) 


for all z, z' £ M d+1 . 


Remark 4.1. An alternative would be to define £ on the sector V = {t £ T: — 2 < 
|r| s < 2 } by the expression 


( n t,X £T )(^ X ) 


where 


[ Q(t — s)(H £ x t)(s,x) ds if |r| s < 0 , 

l T ' (4 - 8) 

[ Q{t-s)(U £ t x T)(s,x)ds + {fl x ,J Q T) if 0 < |r| s < 2 , 

Jt-2T 

(ft,x,J Q T) = ~[ Q(t-s)(U £ tx T)(s,x)ds . (4.9) 

Jt-’2T 


The new symbol is needed to ensure the property (nf x £r)(t, x) = 0 when |r| s ^ 0, 
which is necessary when lifting the fixed-point equation. While this defines an extended 
model with the required properties, the fact that Q acts by convolution in time only limits 
the regularity of its lift to the regularity structure (cf. Remark [477] below). We will see in 
the next subsection why it is sufficient to introduce new symbols £t only when |r| s < 0 . <0 
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Remark 4.2. In cases where v(t,x) takes values in R n and Q(t) G K n , we should in fact 
introduce n commuting symbols £],..., £ n , and define the extended model by the relations 


( n L^d)(^) = 


I 


t—2T 


Qi(t - s)(nLr)(s,x)ds 


(4.10) 


for i = 1 ,... n. The results that follow remain true when £ and Q are replaced by £ ? ; and 

Qi. 0 

We now have to check that (II £ , T £ ) indeed dehnes a model on W. The following lemma 
contains a technical estimate preparing the proof of that fact. 

Lemma 4.3. Assume that for any compact set J4 C M rf+1 there exists a constant C& such 
that 

|(ir j T)(S)|<C*||^-f||l T|l ||r|| (4.11) 

for all t € W and all z,z G A. Then there exists a constant Cq, depending only on Q, 
such that 

\(UI£t)(z)\^C 0 Cj\\z-z\\ 1 s tU \\t\\ (4.12) 

holds for all r G W and all z, z G A, where A = {(t, x): 3t G R: \t — t\ ^ 2 T, (t. x) G J4} 
(for brevity we shall call 14 the 2T-fattening of A, although strictly speaking it is only a 
fattening in the time direction). 

Proof: Using the definition (14.61) of U £ x £t and the assumption (14.111) . we obtain 


l( n t, a A)(*>®)l [_ \Q(t-s)\\\{t,x) - (s,x) 

Jt-2T 


' t—2T 
r2T 


Is ds ||r| 


sc 


r2T r2T 

C-g \Q(s)\\t — t + s\ a / 2 ds \\t\\ + Ca |<5(s)| ds \x — x| a ||r 

Jo Jo 


where a = |r| s . The required bound thus follows if we can show that 

r 2T 


f 

/ |t — t + s\ a ! 2 ds <\t-t 

Jo 


a/2 


holds for \t — 1 1 < 1. By treating separately the cases t > t and t < t, one sees that the 
left-hand side is always bounded above by a constant times (|f — t \ + 2 j i ) 1 + q / 2 ; which is 
bounded above for a > — 2 . Since on the other hand, the right-hand side is bounded below 
by a positive constant for a < 0, the result follows. □ 

Note that m Prop. 8.27] shows in particular that the assumption (|4.11|) is satisfied 
by any canonical model for mollified noise IP built as in Section T3.21 

We can now state the main result of this subsection, which is an adaptation of the 
extension theorem m Thm. 5.14] and of [19] Prop. 8.27] to our degenerate situation. 

Proposition 4.4 (Extension theorem for £). Let Z e = (IF,r £ ) be a model for the reg¬ 
ularity structure (A,W,Q), where W C V. and such that z II £ t(z) is continuous and 
satisfies (14.111) for any e > 0. Let W = W U {£t: t G IP). Then Z £ = (II £ ,r £ ) obtained 
by extending Z £ in the above way is a model for (A,W,Q), which satisfies (14.121) for any 
e > 0. Furthermore, 

l|r e || 7 ; s= inu (4.13) 

holds for any 7 G R and any compact A. C M rf+1 . 
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Proof: We have to prove that Z £ satisfies the assumptions (M1)-(M4) of Definition 13.41 
Property (Ml) is automatically satisfied owing to (14.71) . Regarding property (M2), it 
is known m Sec. 8.3] that there exists a linear map IF : T —» ^(M 0 ^ 1 ) such that 
IF r = II %(F Z ) 1 t is independent of z for any r € T. Property (M2) holds if we can find 
for any r € W a distribution IF(£r), independent of (t,x), such that 

n I x (£t) = U £ F t £ x (£r) Vz = (t, x) € R d+1 . 

To achieve this, we simply define IF on the extended structure by 

(II e £r)(t,x)= f Q(t — s)(ITt)(s, x) ds . 

Jt-2T 

Indeed, we have (writing as usual At = t^ 1 ) <g> t^) 


(n £ F £ z £r)(i,x) 


U e (£® f £ z )Ar(t,x) 

(n'£r (1) )(t,i){/J,T <2, > 


f Q(t - s)(n'r (1 >)(s,i)ds(/J,T |2) ) 

Jt-2T 

[ Q(t - s)(U e z T)(s,x)ds 
Jt-2T 


(n £ z £r)(t,x) . 


To obtain the fourth line, we have used the fact that H £ z = Yl £ F £ holds for the original 
model, and thus 

(n*r)(s,x) = (: U. e z FIt)(s,x) = (II|(Id <g>f z )Ar)(s, x) = (IFt (1) )(s, x)(/f, r {2) ) . 

Property (M3) is a direct consequence of (14.121) . Indeed, this bound implies that for any 
localised scaled test function rj of integral 1 , 


\(Ul x £T,r,l x )\ = 


t — t x — X 


dtdx 

fJ+ ^ 2 s,x + 5y)\\r](s,y)\ ds dy 


<<5 |T| ir|| 

as required by (13.201) . Finally, in order to prove Property (M4), we recall from (14.51) that 


r %£ T = £T lz r • 


If |r| s ^ 0, there is nothing to prove. If /? < a = |t| s < 0, then 

Wn-zZrb = \m- z rb = \\ T lz T \\p < l|r e || 7; s||T||||z “ ^ll?" /3 
for all This completes the proof that Z £ is a model, and also proves (|4.13l) . □ 
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4.2 Multilevel Schauder estimates 

We now would like to construct an operator : £> 7 —>• D 7 which lifts the operation of 
integration against Q to the regularity structure. Formally, this means that we should 
have 

1ZJC® f = K v * 1Zf , K v (t , x , t, x) = Q{t — t)S(x — x) (4-14) 

for all / G £> 7 . In fact, this should really be interpreted as 

(7 lJC®f){t,x)=[ Q(t-t)(Kf)(t,x)dt , (4.15) 

Jt-2T 


or in terms of test functions 

where ip is defined by 




Q(s -t)i/}(s,y) d s 



t)ip{s,y)ds . 


(4.16) 


(4.17) 


The problem with such a plan is that the kernel K v being singular in space, we cannot 
apply an expansion as in (13.40[) and (13.41 [1 . However, for the class of equations we are 
interested in, we do not actually need to define KSj on all of £> 7 . It will be quite sufficient 
to define it on a subset of D 7 , which is given by the functions whose components with 
negative homogeneity do not depend on (t,x). This motivates the following definition. 

Definition 4.5. Let Z = (n,T) be a model. The space T>^(Z) = T> q(T) is the space of 
functions f G D 7 (r) of the form 


f(z)= CtT+ X] c T ( z ) t ='■ f-+ f+( z ) , (4.18) 

rgT : — 2 <| t | s <0 rGT: |r| s ^0 


where the c T do not depend on z. 

Remark 4.6. Not all constant functions belong to as they still have to satisfy the 
analytical bound (13.341) . which includes the requirement \\f(z)—T Z zf(z)\\p = 0(\\z— ^|| 7 _ 4 ) 
as z —> z for all /3 < 7 . In fact, since T zZ t = (Id ( 8 > 7 zS )Ar, we see that a sufficient condition 
for having T zz t = r is 

Ar = r (8) 1 . (4.19) 

Thus the sum defining /_ should involve only terms c T r such that At = r® 1 . As can be 
seen in Tabled! in our case this property is satisfied by the symbols E, T, V, ^, and will 
also hold e.g. for ?, but it does not hold for & for instance. Note that owing to (13.151) 
and (13.331) . (14.191) implies n^r = Ii z r for all z,z G M rf+1 , i.e., the model does not depend 
on the base point z for these r. 0 

With this notation in place, we define the operator on Vq by 

(K.Qf)(z) = £f-+ [ Q{t - s)f + (s,x)ds . (4.20) 

Jt-2T 

Note that this definition indeed requires £ to be defined only on symbols t of strictly 
negative homogeneity. 
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Remark 4.7. An alternative definition of /C 7 , closer in spirit to (13.391) . would be to set 
(£?/)(*) = £f{z) + J Q (z)U(z) + (M?f+)(z) , 
where £ is defined for positive-homogeneous terms as in Remark 14. II and 

n-fc 

J Q (t,x)r= / Q(t - s)(n^T)(s,a;)dsl , 

Jt-1T 

(J\f®f + )(z) = [ Q(t - s)(Kf+ - U^f + (t,x))(s,x)dsl . 

J t—2T 


While it is easily checked that /C 7 defined in this way indeed satisfies (14.151) . the problem 
is that it typically does not map X > 7 into itself, but only into some D 7 with < 1. This is 
due to the fact that we cannot Taylor-expand the kernel K v in space, preventing us from 
subtracting higher-order polynomial terms as in (13.411) . 

The following lemma on translation invariance of the canonical model Z £ = (n £ ,T £ ) 
will allow us to prove that /C 7 given by (14.201) satisfies the required properties. We give 
its proof in Appendix 0 

Lemma 4.8. Let (n £ ,T E ) be the canonical model defined by (13.221) . (|3.27l) and (14.61) . For 
all h, z,z £ M rf+1 , all r £ T, all f £ F + and all e > 0, one has 


K+h r {z + h) = n %t(z) , 

1 z+h,z+h‘ ~ L zz‘ ) 

( 7 l+h,z+h, f ) = ■ ( 4 - 21 ) 

The following result shows that the multilevel Schauder estimates contained in im 
Thm. 5.12] also hold in the case of £, in a similar form. 

Proposition 4.9 (Multilevel Schauder estimates on /C 7 ). Ifj > 0, then the operator 1C® 
maps Vq into Vq, and satisfies 

|||/C?/||| 7; ^(1V||Q|| l1 )|||/||| 7;1 , (4.22) 

f f 2T — 

where IIQH^i = / |<3(t)|dt = / |Q(i)|dt, and 8. is the 2T-fattening of 8. The identity 

Jr Jo 

(TZKSf)(t,x) = f Q(t — s)TZf(s,x) ds (4.23) 

Jt-2T 

holds for all f £ Tfif Furthermore, if (II, T) is a second translation-invariant model 
satisfying (13.271) and (|4.6I) . and /C 7 is the associated lift, then 

|||/C?/ ;7C?/||| 7; * < (1 V \\Q\\ L i)i\f ; /| 7;1 (4.24) 

holds for all f £ 7>q(T) and f £ Pq(T). 
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Proof: We start by showing that relation (14.231) holds, assuming /, f £ Pq. It follows 
from the reconstruction theorem and its corollary |19} Prop. 3.28] (see also ([4. lip ) that 

Kf+(z) = (1 ,/+(*)} = ( IL E J + (z))(z ) 

for all z £ R d+1 . Furthermore, let 

(£/-)(*) := n lf_{z) = Y. 

— 2 <| r | s <0 

(where the right-hand side actually does not depend on z, cf. Remark 14.61) . Then we 
see that TZf- trivially satisfies the reconstruction theorem, so that by uniqueness of the 
reconstruction operator if 7 > 0 we have 

Kf- = n.f- = n s j_ 

whenever /_ £ Pq. Thus in fact 

Kf(z) = {Ulf(z))(z) 

holds for all 2 : £ R rf+1 . It follows by (14. 5p and the same argument as the one yielding 
TZf- = FI 2 /_ that 


(K£f-)(t,x) = (n l x £f_)(t,x) 

= I Q(t~ s){Ul x f-)(s,x)ds 

J t—2T 

= f Q(t — s)lZf-(s, x) ds , 

Jt-1T 

which implies that we indeed have 

(7 = f Q{t-s)Kf + (s,x)dx 

Jt—2T 

= j Q(t — s)1Zf(s, x) dx . 

Jt—2T 


Next we check that /C 2 indeed maps Pq into itself. Since £ preserves homogeneity, it is 
in fact sufficient to show that maps Pq into P 7 . First we note that for any (3 < 7 , 

\\K$f(z)b^W-(z)b+ I IQ(i-'S)|||/+(s,x)|| /3 ds 

Jt-2T 

< \\f-( z )b + IIQIU 1 su p \\f+( s , x )b 

se[t-2T,t] 

^ (1 v ||<5 I|l i )III/III 7; i > 

(note that depending on the sign of /3, either one or the other term in the sums contributes, 
but never both). Regarding terms evaluated in different locations, note that (14.51) implies 
that 


11 si- - n- z £f-b = \ m ~ - n s f-)h 

= \\f--r%zf-b 

< i/-iuk - *nr* 
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As for the positive-homogeneous part g = /C 7 /+, using Lemma fTTBl we get 


g(t, x ) - r E tiX .i tS g(t, x ) 


and thus 


r2T 

/ Q(«)[/+(t 

J 0 

i-2T 

/ Q(s)[/+(* 

J 0 


a >*) - r t,x;f,s/+(* -*>*)] ds 

»>*) - a >®)] ds » 


r2T 

II 9(t,x)-Tl is g(t,x)\\p < / |Q(s)|||/+(i-s,a:)-rf f / + (t-s,x)|| /3 ds 

Jo 

Collecting the estimates obtained and using the definition (13.341) of jj| • |||, we see that 
/C 7 / € D 7 and that it satishes (|4.22f) . 

Consider finally the case where we have a second function / E Dq (U)• Using again (14.51) , 
it is immediate to see that 


|£/- ; ^ III/- 5 /-ll 7; I • 

Writing g(z) = K.® f + (z), it follows in the same way as above that 

\\9-9\\ r ,& ^ ll ( 3lk 1 ll/+ - /+II 7 ,I • 

Finally, we have 

g(t, x) - g(t , x) - r^. t - s 5 f(i, x) + T t}X .^g(t, x) 

/»2T 

= / <9(s) [/+(* -s,x)~ f+(t - s, x) - rf - s, x) + r - A *)] ds • 

Jo 

Using translation invariance of r e and T, we obtain for all j3 < 7 
IKiu) - <j(*>&) - rf )a:;f)2 5(t,x) + T tiX .i tS g(t,x)\\p 

< IIQIU 1 SU P ||(t-s,x) - (t-s,x)||3 ,_ ^|/ + ;/ + || ^ 

se[0,2T] 

= llellL.|l((,x)-(f.x)||r' , |/ + ;/ + || T 3i, 

from which (I4.24|) follows easily. □ 

4.3 Extension to the space D 1,ri 

As already mentioned at the end of Section 13.31 the definition of the spaces (D 7 has to be 
modified in order to be able to deal with the time-zero singularity of the heat kernel. In 
particular, we have to check that the estimates on /C 7 established in Proposition 14.91 still 
hold in these modified spaces. 

We recall a few notations from [19], Section 6 ]. Let P = {(£, x) G M d+1 : t = 0} denote 
the time-zero hyperplane, and introduce the scaled norms 

||z||p = 1 A |f|U 2 , || 2 , z||p = ||z||p A p||p = 1 A |f| 1//2 A \t | 1//2 . (4.25) 

For a compact .h C M d+1 , let 

Up = {(z,z) 6 (R\P) 2 : z^z, \\z - z\\ s ^ ||z,z||p} . (4.26) 
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Definition 4.10 ( [191 Def. 6.2]). Given a model (II, T) and constants 7 > 0 and r/ £ M, 
set 


7 ,??;^ := sup sup 


WMWp 


ze&\p /3<7 ||2;||p ^ )A0 ’ 


n^r 11 / 0 ) 11/3 

Uh,V;& : = SU P SU P I, | ,n_fl 
zeR\P / 3<7 Ill'll p 


(4.27) 


The space V' y,ri consists of all functions f : M d+1 \ P —> T such that 


7 , 7 /;^ 


7 ,wF 


+ sup sup 


\\f(z)-T z -J(z) 


= 117-^1 


Z, Z 


-IIV-7 


< OO 


(z,z)efip /3<7 \\z — z\\ s 
Given a second model (II, T) and f £ V^^lT), we also set 

\\f(z)-f(z)-V zz f(z)+T zz f(z) 


(4.28) 


III / ; / 1 W ;S := 11/ - /|| 7 ,*j;* + sup sup 

(z,z)e&p j3<i 


z — z 


= 117-^1 




^tl’7—7 


(4.29) 


Finally, similarly to Definition |^.<5[ we define Dq ,7? to be the set of f £ T> 7,r? whose com¬ 
ponents with negative homogeneity do not depend on z. 

Most of the results in [19) Section 6 ] are directly applicable to the present setting. The 
only result that has to be adapted is Proposition 6.16, which takes here the following form. 

Proposition 4.11 (Multilevel Schauder estimates on Dq’ ?7 )- Let f £ V^iY 6 ) with 

-2 <r/<0 <7 <r/ + 2, (4.30) 

and assume that f+(t,x ) = 0 whenever t < 0. Then K.®f £ T>q , 7] (Y e ), and there exists a 
constant C\, depending only on Q, such that 


[1 + ||r'|| 7 .j,] l/L,.*, 


(4.31) 


where & is the 2T-fattening o/H. Furthermore, let f £ 'D^'fY), where (II, T) is another 
translation-invariant model satisfying (13.2711 and (14.61) . and assume that f+(t,x ) = 0 
whenever t < 0. Then there exists a constant C 2 , depending only on Q, such that 


|||/C?/;^/||| 7 ^ < Cr[l + l|r £ || 7 ; 1 ]|/;/ 1| 7iT?;1 + C 2 ||/|| ml ||r-r £ || 7 ^ . (4.32) 

PROOF: It follows as in the proof of Proposition 14.91 that 

\t£f-^,w,z < I/-I™* 7 |£/_ < III/- ;/-lll 7 ,„;i, 

so that we only have to consider the positive-homogeneous part /+ of /. For convenience, 
we again write g(z) instead of (/C 7 f + )(z). Note that 


\\f+\\y,v;& = sup sup 


\\U(z)b 

z£R\PP<7 H^llp ^ A ° 


= sup 


sup 


z=(t,x)£$\P 07/3<7 


(1A |t|)(' 5 -’>>/ 2 II/+MII/J 


since we may always assume (3 ^ 0 and thus (7 — /3) A 0 = 7 — /3. Since f+(t, x) = 0 for 
t < 0 we have 

\\9{t,x)\\/3 ^ [ \Q(t-a)\\\f + (s,x)\\pd8^ [ ds ll/+ll 7 , 7 ?; g • 

J (t—2T)V0 JO (1 A S)(P 
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Since /3 < 7 < r/ + 2 , the integral is convergent and we have 


where Cq depends only on Q. In order to estimate |||g||| 7)??; £, we also need to bound 
differences of the form g(z) — T £ zz g(z). From now on we assume that z = ( t,x ) and 
z = ( t , x ) with t ^ t , since the case t > t follows by symmetry. Using the fact that f + is 
supported on {t ^ 0 } and translation invariance, we get 

/*2T At 

II g(z) - T lz9(z)\\p ^ / \Q(s)\\\f+(t - s,x) -T% { s f + (t-8,x)\\pds . 

J 0 

However, there is a difficulty because (z, z) £ &p does not automatically imply that one 
has ((t — s, x ), (t — s, x)) £ £p. In fact, 

((t — s,x), (t — s,x)) £ &p (1 A \t — s| A 1 1 — si ) 1 / 2 ^ ||(i — s,x) — (t — s,x)|| s 

1 A \t — s| A \t — s| ^ \\z — z\\ 2 . 

By definition of Sip, we can assume that ||z — z\\ s ^ 1, so that s ^ t — \\z — z\\ 2 is a 
sufficient condition for having ((t — s,x), (t — s,x)) £ 8.p. We split the integration interval 
into [0, t — \\z — z\\ 2 } and [t — \\z — z || 2 , 2 T At], and treat each interval separately. For the 
first interval, we use the definition of |||/+||| to get 


rt—Wz—zW 


IQ(s)lll/+(* ~s,x)~ r ~ s,®)ll^ds 


rt—\\z—z\\ 


< 




__fls II r _ ^||7-/3||| f III 

(lA(f-s))H )/ 2 11 lls lU+ h,n;& 


c t 


Q 


(1 A t )( 7-»?)/2 


\z-z\\r p 11 


117 , 77 ;^ * 


For the second part of the integral, we bound ||/+(t — s,x) — Y £ zz f + {t — s, x)\\p above by 
the sum ||/+(t — s,x )|| / 3 + \\T E zz f + (t — s,x)\\p, and treat each term separately. Using the 
assumption f+(t,x ) = 0 for t < 0 and performing the change of variables s = t — s, the 
first term becomes 

|Q( t - 5 -)|||/ + (»,x)|| S d»< jf (1 7 -U ll/ + ll 7 „;I 

< ik - iiz+ii,.,,! • 

For the second term, we use the definition of ||r £ || and the decomposition of /+ as a linear 
combination of r to obtain 


|r*,/ + (s,i)l|„ « ^l|/ + ( 5 -,x)||j|k - 2||Ul|r'll T « 

5 

<y Wf+hrr$ || -ll^llrui - 

"" ^ (! A s )( 5_r ?)/ 2 Z lls 11 
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where the sum runs over finitely many values of 5 £ [0, 7 ]. Setting s = t — s, this yields a 
contribution 


l 


\z-z\\i+t-t 


I Q(t - s)|||r^/ + (s,x)|| /3 ds 

^E 


m-s)\ 

(lAl)(H/2 as||Z " lls 


11/4 


-iir e n - 


< 


El4- 2 1. 2_i+ "l4-^li;-' J ll/+ll,,,;5llv 


1 7;J? 


< \\ z _ z\\ 2+r i-P 
^ II^ ^lls 


-iir £ ii - 
117 ;^ 


Combining the different bounds, we obtain (|4.31l) . 

It remains to obtain a similar bound on I#; g |||, where g = /C^/ + . The same argument 
as before yields 

I Is - 9\U,mR < CqII/+ - /+ll 7) „ ; j* • 

Regarding the second term in the definition of |||g ;^|||, we again split the integral defining 
it into two parts. The first one is 


L 


t- z —2 


IQ( s )lll/+(* -s,x)~ f+(t -s,x)~ r %f+(t - s,x) + T z ~f + (t - s, x)\\p ds 

Cq 


St 


(1 A t)(T'-^)/2 


^-2|ir' J l/+;/+ll 




For the second part, we use the decomposition 


||/+(f -s,x)~ f+(t -s,x)~ T e z J + {t- s,x) + T zS f + (t- s,x)\\p 
||/+(* ~s,x)~ f+(t - s, a:)||/3 
+ IIn 5 (/+(t- s,x) - f+(t- s,x))||^ 

+ \\(Tl- z -T zz )f + (i-s,x)y, 


and estimate the integral of each term separately, in the same way as above. Combining 
the different bounds, we obtain (14.321) . □ 

Remark 4.12. Note that the condition (14.301) on 7 and 7 excludes integer values for these 
parameters, except that we may have either 77 = —1 or 7 = 1 . The reason while these 
integer values are allowed is that owing to the parabolic scaling, the exponents such as 
(7 — 77) /2 occurring in integrals over time can never take integer values that could make 
them divergent. 


5 The fixed-point equation 

In this section, we construct a fixed-point equation lifting (14.21) to an appropriate space 
X> 7,7? , and prove that it admits a unique fixed point, depending Lipschitz-continuously on 
the lifted initial data and noise. The procedure is very similar to the one for the Allen- 
Cahn equation or the 4>3 model carried out in m Section 9], with only a few adjustments 
required by the additional variable v. 
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5.1 The set-up 

From now on, we assume that the space dimension is d = 3, since this is the most difficult 
case which is still renormalisable. Consider the fixed-point equation on T)^ ,ri x given 
by 

U = (lCy + R^U)R + [H + F{U, V )] + Gu 0 , 

V = JC®R + U + Qv 0 , (5.1) 

where R + (t,x) = l{t>o} and 7 ^ 7 — 2 as before, and by definition, IZR+z, is given by 
the distribution £ e l{t>o}- Furthermore, 

(Gu 0 )(t,x) = (, S(t)u 0 )(x ) = [ G(t, x -y)u 0 (y)dy , 

J R 3 

(Qvo)(t,x) = e tM «o(a;)l{t>o} • (5-2) 

In fact, these distributions have to be lifted to D 7,7? , which is always possible by con¬ 
struction, cf [T9) Def. 2.14], In order to understand the structure of the solution of the 
fixed-point equation (15.11) . let us rewrite it in the form 

U = 1R + [E + F(U, V)] + R(U, V ) + Gu 0 , 

V = £R+U- + QR+U + + Qv 0 , (5.3) 

where U- and U+ denote respectively the strictly negative homogeneous and positive 
homogeneous parts of U (cf. (|4.18|> ). R(U,V) contains the polynomial part of /C 7 [E + 
F(U, 14)] as well as the term i? 7 72. [E + F(U , V)] (which always belongs to T), and 

Qf+{t,x)=[ Q(t - s)f+(s,x)ds . (5.4) 

J t—2T 

Assume F is a cubic polynomial of the form 

F(u, v ) = a\u + a 2 v + /3i u 2 + (3 2 uv + P 3 V 2 + 71 u 3 + r y 2 u 2 v + 73 uv 2 + 74 U 3 . (5.5) 

Further assume that the initial conditions rto and vq are such that Guq and Qv 0 are 
sufficiently regular. Iterating the map (15.31) . starting for instance with identically zero 

functions, it is not difficult to see that any fixed point of (15.11) must have the form 

U = T + cpl + [ai Y + a 2 Y + a 3 Y + a 4 Y] + [&i Y + b 2 Y + 63 Y] + (V<p,X) + qu , 

V = ? + -01 + [ai Y + a 2 Y + a 3 Y + a 4 Y] + [& 4 Y + b 2 Y + 63 Y] + (V^, X) + g v , 

(5.6) 

where the coefficients satisfy Oj = 7 j for i = 1,2, 3,4, 

bi = Pi + 3<£7i + ^72 , 

b 2 = fi 2 + 2(^7 2 + 2^73 , (5.7) 

bo = h + ^73 + 3^74 , 

and diT = Q(ajr), 6 jT = Q( 6 jr). In (15.61) . gu,Qv denote terms of homogeneity at least 
~ — Ac, and (/?(z), ^( 2 ), Vv?(z), Vt/>(z) are independent functions — the notation is not 
supposed to suggest that f and ip are differentiable. The scalar product notation (Vc p, X) 
is a shorthand for Yli=o(^F)i^i- 
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Proposition 5.1 (Equivalence of fixed-point equations). Assume thate > 0 and that (15.11) 
admits a fixed point (U, V ) where U, V £ D 7,7 for some 7 , ?? £ M. Then ( u , n) = (JZU, TZV) 
satisfies the fixed-point equation AMD- 

Proof: The first step is to note that (13.271) and (14.61) imply 

( n ! T ) 0 ) = J K i z ~ 2 i)£ £ (zi)dzi =: Xe(z) , 

( n l ? )(^) = J I<Q ( Z - zi)i e {zi)dzi =: Xe( z ) > 

independently of z, where we have used the notation 

/»2T 

K®(t,x)= / K(t — s, x)Q(s) ds . 

Jo 

Note that both in dimensions d = 3 and d = 2, the sum in (13.271) is empty, because in 
both cases T has strictly negative homogeneity, cf. Table [H Applying the reconstruction 
operator 1Z to (|5.6D . since all functions are smooth for e > 0, we obtain 

u(z) = Xe{ z ) + T( z ) , 

v( z ) = X?(z) + ip( z ) ■ 

Here we have used the fact that if /+ £ has only components of strictly positive 
homogeneity, one has 7Zf + (z) = (U £ z f + (z))(z) = 0 (cf. (|4.11j) as well as (T9], Prop. 3.28]). 
On the other hand, using (|3.38l) and (13.431) . we obtain for all / £ P 7,7? 

{Kite + RyK)f) (z) = (G * TZf)(z) . 

Applying 1Z to the equations (15.11) and using (|4.23l) . we thus obtain 

u(t, x) = G * 1ZR+ [H + F(U, V )] (i, x) + Gu 0 {t, x) , 

v(t, x) = f (7ZR + U)(s,x)Q(t — s) ds + Qvo(t,x) . 

Jo 

It is not difficult to check that 1Z and R+ commute. Thus the result follows if we are able 
to show that 

1ZF(U, V ) = F(1ZU, 1ZV) . 

Computing the part of nonpositive homogeneity of the difference 1ZF(U, V ) = F(1ZU, TZV), 
we see that it is a linear combination of terms of the form 

(H^)(z), (n^)(z) and (ffiT)(.~) 2 (V^, (H;I)(z)) , 

and all possible terms obtained from these by substitutions of the form T >—)■?. It follows 
directly from (|3.22l) that (H;A)(z) = 0. It thus remains to show that all terms of the 
three other types vanish as well. Applying A3.271) . A3. 281) . and recalling A3.14|) we obtain 

(n§ Y)(*) = J[K(z - Zl ) - K(z - Zl )] (n§ V)(*!) dz !, 

which vanishes in z = z. It follows that (n|*Y*)(z) = {U. £ z Y) (z) (n| V) (z) = 0. All other- 
terms can be treated similarly, and the result follows. □ 
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Remark 5.2. In dimension d = 2, the symbols have strictly positive homo¬ 

geneity, so that the last part of the proof is not needed. 

The expansion (15.61) shows in particular that both U and V' have regularity — ^ — n, 
and that their components of negative homogeneity do not depend on (t,x). We will thus 
look for a fixed point of (15.61) in a space V g ,T? for appropriate values of 7 and 7 . 


5.2 Local existence and uniqueness of the fixed point 

In order to describe the effect of F, it will be useful to introduce the following slight 
generalisation of Definition 14.51 Given a model (II, T), the space Vj ,r ^(T) is the space of 
functions / € V' y,ri (T) of the form 

f(z)= ^2 c t t+ ^2 c T (z)r , (5.8) 

t€T: a7Ms</3 tGT: |r| s 7/3 

where the c T do not depend on z. 

The following result establishes the strong local Lipschitz continuity of F in the sense 
of [19], Section 7.3]. 

Proposition 5.3 (Strong local Lipschitz continuity). Fix a < 0 < 7 and set 7 = 7 + 2 a 
and fj = 3rj A (7 + 2a). Let It C M d+1 be compact, and assume U, V € Pq’^(T) are such 
that 

( 5 - 9 ) 

Then F(U,V) eV™ (T), and there exists a constant C(R) such that 

I F(U,V)\l^^C(R) . (5.10) 

Furthermore, if (II, T) is a possibly different model and l/,f £ Dg^(r) satisfy the bound 

III^IIIt, v;& III ^lll'y,n;^ ^ 


III F(U, V ); F(U, V) lly,*,;# ^ C(R) |||U ; U\\^ + 


117,»? 


■,& + ll r ~ •^’ 1127 +, 


0+ 


(5.11) 


Proof: Writing U(z) = U- + JJ+(z) where U- contains all terms of negative homogeneity 
and applying m Prop. 6.12], we see that 

U(z) 2 = U 2 + 2 U-U+(z) + U + {zf € V l+ a a 2 ^{v+o<) ^ 

U(zf = U 3 _ + 3 U 2 U+{z) + 3 U-U(z) 2 + XJ + (zf G V^’ 3vMr,+2a) , 

where the multiplication of elements in T>?f is well-defined due to the results in [ 191 Sec. 4]. 
Similar relations hold for powers of V and cross-terms, proving (15.101) . The proof of (15.111) 
is similar, using the bound on |||/; (/||| given in [19| Prop. 6.12], □ 

Remark 5.4. In cases where the two models (II, T) and (II, T) are identical, (15.111) auto¬ 
matically provides the bound 


Him v) - F(U, V)h m < C(R) III U - U ||| 7 ,^ + |||P - V\\^ 


(5.12) 
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In order to proceed, it is convenient to rewrite the fixed-point equation m in a 
slightly different form. First we set 

W := (/Gy + R^R+E . (5.13) 

For the time being, let us assume that W £ ^o’ao +2 f°r some 7 G R and all rj < Qo + 2 
(recall that ao = — | — k in dimension d = 3); we come back to this point in Remark 15.61 
below. Then we define the map 

M : U ^ (/Gy + R 7 ^)fi+F(C/, JC^R+U + Qv 0 ) + IF + Gu 0 . (5.14) 

If U is a fixed point of At and V = /C 7 R+U + Quo then ([/, V) is a fixed point of (15.11) . 

Following [T9l Section 7.1], we write O = [—1,2] x R 3 and Ot = (—oo,T] x R 3 , and 
introduce the shorthand ||| • | 7 ,r/ ; T for |j| • || 7 jT? ;O t , and define ||j • ; • | 7 , 7 ?; r similarly. Finally, we 
use the notation \Z ; Z||| 7; ^ = ||II — n|| 7; ^ + ||F — r|| 7; ^ to quantify the difference between 
two models. The following result is an adaptation of m Thru. 7.8]. 

Proposition 5.5 (Existence and uniqueness of the fixed point). Assume — | < r? < — ^ 
and a > — 1 are such that r] + 2a > —2, and assume 7 > —2a. For any uq,vq,W such 
that IF + Guq,Qvq £ V™, there exists a time T > 0 such that At admits a unique fixed 
point U* £ V y,ri on (0, T). Furthermore, the solution map St ■ (uo,vo,W, Z) i->- U* is 
jointly Lipschitz continuous in the sense that if Z and Z are two models and (uo,vo,W) 
and (uo,vq,W) are two sets of initial conditions and forcing terms such that 

|||Guo; Guq ||| 7 j 77;t + HI+ |||VF; JF||| 7!7?; t + \Z ; Z|2 7 + a; o ^ <5 > (5.15) 

then one has 

|l/*;Cr*ffl 7 j „;T<C , o(J (5.16) 

for some constant Co > 0 . 

Proof: Given R > 0 let B(R) = {/ £ Dq’^: ^ R} and pick U £ B(R). By 

Proposition 14. 1 ll we have 

V := K%U + Qv 0 £ B(C[R + |Qu 0 | 7 iW t) 

for a constant C[ depending only on Q and on the model. Applying Proposition 15.31 we 
see that F(U, V ) £ R^oiza an d satisfies 

\\F(U, V)\l^, T st C{R + C[R + IQvoIwt) • 

The assumptions on 7 , q and a guarantee that 7 > 0 and fj > —2. We can thus apply 
Theorem 7.1 and Lemma 7.3 in m to obtain the existence of constants C 2 ,k > 0 such 
that 


111^(^)111^ C 2 T*/ 2 C{R + C[R + |Q«o| 7 > „ ; t) + II^| 7 , 7?; t + \lGuo\U, V ;T 

for every T £ (0,1]. Choosing 

R = i+OI^IUt + IIIGuoIIUt] 

and T sufficiently small, we see that At maps the set B(R) into itself. 
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Next we show that M. is a contraction in B(R). For U, U £ B(R ), V as above and 
V = /C 7 Z7 + Qu 0 , we have 

M{U)-M(U) = {IC^ + R^n)R+[F(U,V) - F(U,V)\ . 

Using again Theorem 7.1 and Lemma 7.3 in |19| and applying (15.12[) to bound the term 
[F(U,V)-F(U,V)\, we get 

\\M(U) - M(U)%, mT < C 2 T«' 2 C{R + C[R + |||Quo||| 7 )??; t) [||U - V\ lt ^] . 

Taking again T small enough, we obtain that Ad is indeed a contraction, and the existence 
of a unique fixed point in a ball in T > 7,,? follows from Banach’s fixed-point theorem. 

Finally, in order to prove the bound (15.161) . we set U\ = (/C 7 + R~f1Z)R + F(U,V) and 
U i = (/Gy + R~/7l)R + F(U , V). Using the inequality 

1/ + 5; / + < 1/; + Ills 1 ; 

we obtain 

\\M{U)-M(U)\\^,t ^ |||Ui; UiII^t + \\W+ IGuo-,Guoh, m T , 

where the first term on the right-hand side can be estimated with the help of (15.111) . The 
required bound thus follows from pT9l Thm. 7.1 and Lemma 7.3]. □ 

Remark 5.6. In fact, Proposition 15.51 can be applied as soon as we choose 7 > 1. Indeed, 
let a = — \ — k and assume 7 = 1 + R for some R > 0. 

1. The conditions on 7 and 7 are satisfied provided k < R and k ^ i. 

2. Lemma 7.5 in [19] shows that if uq £ C v (with the Euclidean scaling and recalling that 
we take —| < 7 < —|) then Guq € P 7,r? for any 7 > (7 V 0). 

3. Definition 2.14 in m shows that if vq £ C 7 , then it can be lifted to a an element 
uo £ D 7 , and then Qv 0 will belong to X >7,J? as required. 

4. As discussed in Section 9.4 of [19] (see Equation (9.18) and below), the assumption 
W £ P ™ +2 is indeed satisfied for 7,7 as in Proposition 15.51 provided we define 
1 ZR+E as the distribution £ £ 1 {t^o}- 0 

Let Z £ = T(£ £ ) denote the canonical model for mollified noise £ £ constructed in Sec¬ 
tions 13.21 and ED In exactly the same way as in m Prop. 9.8 and 9.10], we can now 
extend the solution map St to a maximal solution map S L , defined until the first time 
\\(TZU)(t, •)]],, + ||(7 ZV)(t, •)||j ? reaches a (large) cut-off value L. This map satisfies 

HS L (u 0 ,v 0 , Z £ ) = S L (u 0 , u 0 , ) , (5.17) 

where S L is the solution map of the original SPDE-ODE system (14.11) with mollified noise 
£ £ . In other words, the following diagram commutes: 

(u 0 ,u 0 ,Z £ )i SL > (U*,V*) 


S L 


{uo,v 0 ,£ e ) h 
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(u £ ,v £ ) 












6 Renormalisation of the equation 

6.1 The renormalisation group 

Consider a general model Z = (II, T) for our regularity structure, which is admissible in 
the sense that is satisfies (13.271) . the last identity in (13.281) and (14.61) . The model may also 
be specified by the pair (II, /) such that for all z, z € M“ +1 , one has 

n z = nr /2 and T zg = (T/J -1 ^ =: r 7z2 . (6.1) 

To build a renormalisation transformation, one first needs to define sets T* C To C Tf 
such that 

• To contains all r £ Tf of strictly negative homogeneity; 

• for all r € To, one has At £ span(To) <8> span(To~), where Tq~ C Tp contains all sym¬ 
bols of the form (13.111) with Tj £ T* and \J'k j Tj\ s > 0 (including the empty product). 

See [HI Sect. 8.3], in particular Remark 8.37, for a proof that it is always possible to find 
sets J 7 * and To satisfying these two properties. For the Allen-Cahn model, a possible 
choice is 


t; = {^, v,T}, 

To = {3, V, V, T, «S, VJi, Y-, *^,1, 4-, Y, V,Xi} . (6.2) 

For FitzHugh-Nagumo type equations, we choose the same sets, enriched by all symbols 
obtained by substitutions !—>•?. 

Let T~Lq = span(To) and Hq = span(To“). By definition, a renormalisation transfor¬ 
mation is a linear map M : T-Lq —> T~Lq such that MX k = X k , MIt = TMt for all r £ To 
such that It £ To and M£t = £Mt for all t £ To- With M, we would like to associate 
another admissible model (II M ,/ M ) by setting 

U m t = UMt Vt £ T 0 . (6.3) 

As shown in Till Section 8.3], this is indeed possible if one is able to find linear maps 
A A/ : Hq —> Hq <8> "Hq, and M : —y Hq such that 

nf t = (n* <g) f z )A M r , 

(/f,r) = (/ z ,Mf) (6.4) 

holds for all t £ To and all f £ Tq~, where IlY and ff r satisfy relations analogous to (16.11) . 
The model ( U M J M ) is admissible provided the conditions 

MJk = M(J k <8> Id)A M , 

(Id®A4)(A®Id)A M = (M®M)A (6.5) 

hold, where A4 is the multiplication map, defined by = t<j. The map M should 

also be a multiplicative morphism leaving the X k invariant, that is, 

M(tit 2 ) = (Mt!)(Mt 2 ) , M(A fc ) = X k . (6.6) 

It is shown in [HI Prop. 8.36] that given a linear map M, there is a unique choice of maps 
and M such that (16.51) and (|6.6D hold, making (II M ,f M ) an admissible model. In 
addition, the elements Tf| can be built using a map A M : Hq —> T~Lq <S> Tq such that 

^ = {T~ z ®h)A M . (6.7) 
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Definition 6.1 ( {!?, Def. 8.41]). The renormalisation group 91 is the set of linear maps 
M : TLq —> hio satisfying MX = XM, MS = SM, M{X k ) = X k , and such that for all 
t € X ,o and all t € Fff one has 

A m t — t (g> 1 € Q)(Tp n span^o)) ® T + , 

h>a 

A"f - r ®i e ®(r s (~1 span^Q - )) ® T + , (6.8) 

h>a 

for all r £ T a and all t £ T a C I span(J r ( ) h ). 

It is shown in [19, Lemma 8.43] that 91 is indeed a group. The main result [19., 
Thm. 8.44] states that given M € 91, and an admissible model (II, /), the model (IT 17 , f M ) 
constructed as above is indeed an admissible model. In particular, it satisfies the required 
analytic bounds and behaves well under the extension theorem [H9l Thm. 5.14]. 

Let us now introduce a particular subgroup of 91 suited to our problem. In the case 
of the Allen-Cahn equation (cf PE] Sections 9.2 and 10.5]), a two-parameter group of 
transformations M = exp{— C\L\ — C 2 L 2 } is sufficient, where the generators L\ and L 2 
are defined by the substitution rules 

Li : V H 1 , L 2 : ^ ^ 1 . (6.9) 

More precisely, Lit is in general a sum of terms, where each term corresponds to one 
occurrence of V = X(H) 2 in r. Thus for instance Li(^) = 3T, because there are 3 ways 
of picking a pair of T in , and replacing any of these pairs by 1 yields T. 

In our case, we may a priori have to apply similar substitutions for all elements ob¬ 
tained from V and by replacing one or several L(H) by SXifS). We thus introduce 
the sets of symbols 

■Fi = {V, V, V}, 

X 2 = {L(ti)t 2 : n,T 2 € 7i} . (6.10) 

Note that X 2 has 9 elements. This leads us to define a 12-parameter subgroup of 91 given 
by 

M = exp j— ^2 C\(t)L t - ^2 C' 2 (t)L t | , (6.11) 

r£Ji rgT2 

where L r denotes the substitution r 4 1, applied as often as possible. The exponential 
is then defined by its formal series, which is always finite because any application of L r 
strictly decreases the number of occurrences of X(H) or £(Z(H)). In practice, however, we 
will see that only 2 of these parameters are really needed to renormalise the model. 

We introduce the shorthands C\ = C'i(V), C[ = Ci(V) and Cf = Ci(V). The 
action of M on a few elements that will occur in the computations below is given in the 
following list: 


M(T)= T , 

(6.12a) 

Af(V) = V -CU , 

(6.12b) 

M(*V) = 'V — 3C*i T , 

(6.12c) 

M (— C 2 () 1 — Ci *T* , 

(6.12d) 

M (— 3C 2 ( ) T — 3(7i *st* — Ci T + 3 Cf ''' . 

(6.12e) 
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The action of M on elements such as V is obtained by obvious substitutions, e.g. we have 
M(^) = ^ — CfT — 2(7{?. As in [19j, Section 9.2], one can determine the map A M and 
thus, via (16.41) . the renormalised model II M . Again, we just list a few expressions that will 


play a role in the computations below: 

nf (T) = n z (T) , (6.13a) 

nf (v) = n.(T) 2 -Ci , (6.13b) 

nf(^) = n 2 (T) 3 -3Cin 2 (T) , (6.i3c) 

nf (Y) = n z (Y) , (6.i3d) 

nf(-Y-) = n 2 (Y)nf(v)-c 2 (^) , (6.i3e) 


nf(-Y-) = nf(Y)nf(v)-3C' 2 («Y*)nf(T) + 3C' 1 n 2 (vx i )(/„j ei (T)). (6.i3f) 

Similar relations with obvious substitutions hold for the elements obtained by having 6 
act on some Z(H). 

6.2 Convergence of the renormalised models 

From now on, we assume that the driving noise £ is Gaussian white noise on D = M x if 
where T d = is the d-dimensional torus . This means that we are given a probability 

space (11, J 7 , IP), the Hilbert space H = L 2 (D), and a collection of centred jointly Gaussian 
random variables {Wh}heH which are assumed to satisfy 

E {W h W~ h } = {h,h) (6.14) 

for all h,h £ H. Then £ is the distribution defined by (£, ip) = W nip for every test function 
ip, where ir is the canonical projection from L 2 (M rf+1 ) to L?(D). 

Let Z e = (n e ,T e ) be the canonical model built for mollified noise £L Our aim is now 
to find a specific sequence of renormalisation maps M £ in 91 such that the sequence of 
models Z £ = M e Z £ defined by 

nf = (n £ j M * (6.15) 

converges in a suitable sense to a limiting model Z. 

The relevant result to achieve this is Theorem 10.7 in m • It states that if for all 
localised scaled test functions rf z = S d z r] and all r € T- = {r £ Tf'- |t| s } < 0, one can 
find random variables n,r such that 

E\(U z r,r, 5 z )\ 2 <S 2 ^ +K , 

E|(n 2 r-nf r,7/f)| 2 < e 2 VM*+« (6.16) 

holds for some k, 6 >0, then there exists a unique admissible model Z such that for all 
compact sets ^ C M d+1 and all p ^ 1 one has 

E||Y!^<1, Mill Z-Z £ f rJi <E ep . (6.17) 

Here |||Z||| 7; ^ = ||n|| 7; ^ + ||T|| 7; ^, and |||Z ;Z||| 7; ^ is defined in Section lA2l In (16.161) . each 
random variable n 2 r should in addition belong to the inhomogeneous Wiener chaos of 
order ||r||, where ||r|| is the number of occurrences of E in r (cf [T9. Section 10.1]). More 
precisely, Wiener’s isometry allows one to define linear maps Ik '■ H® k —>• L 2 , as explained 
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in m Section 10.1]. Then the Wiener chaos expansion of IIq r is obtained by finding 
distributions (W’^ £;fc V)( 2 ;) G H® k such that for every test function <p, 


(n ( 0 £) r,<p) = 


S Ik 


fe^ll 


v{z)(yv^ k \){z)&z 


(6.18) 


Let further Vy( £ ; fc ) be the distributions defined in a similar way from the bare model ILL 
As usual, we will express (16.1811 by the slightly more informal notation 

(flS'V)W = XI h{(W^T)(z)) . (6.19) 

/c^||r|| 


Proposition 10.11 in [19] provides conditions on the W^ £;fc ) from which the required con¬ 
ditions (16.1611 follow. Namely, for each r G J 7 -, one should find functions such 

that 


((W^t)(z),(W^t)(z)) J2 {\\4* + M*) \\z~z 

Bk,r 

^Ce 26 ^2 (Iklls + ||2||s) A p - 


^11-v ^ ||K+2 |t|s—A 

s 


((5W^t){z),(5W^t)(z)) 


A || ^ _M k-\-2\t\ 5 — A 

z \\s 


( 6 . 20 ) 


for some k, 9 > 0, where SW^ S]k ^ = and each Bk jT is a finite set of indices 

A G [0, 2|t| s + k + |s|), where |s| = d + 2 = 5. Furthermore, the scalar product in (|6.20l) is 
the canonical scalar product on H® k induced by the scalar product on H. 

In order to obtain such bounds for our system, we can to a large extent adapt the 
proof of m Thm. 10.22]. As a matter of fact, it is immediate to see that this theorem 
can be extended to our case, because the involved kernels have singularities which are not 
worse than in the Allen-Cahn case. The aim of the discussion that follows is thus only 
to determine which parameters of the renormalisation group introduced above are really 
needed to ensure convergence. 

The first step is to observe that if g £ is the mollifyer, i.e., if £ £ = g £ * £, noting 
K e = K * g £ we have 


(n§T)(z) = ( K E* 0 ( Z ) = J K £ {z- zi)£(zi)dzi =:Xe{z) , (6.21) 

which is independent of z and belongs to the first Wiener chaos with 

((W^l)(z))( Zl ) = K £ (z - Zl ) . (6.22) 

As a consequence, we simply set 

(n? } T)(*) = Xe(z) , (n*T)(«) = xo(z) := (K*£)(z) . (6.23) 

Due to translation invariance, it will henceforth be sufficient to evaluate models in z = 0. 
As already seen in the proof of Proposition 15.11 we have 

( n o ? )(X) = ( kQ *£ e )00 = J K Q (z- z 1 )^ E {z 1 )dz 1 =:Xe{ z ) ■ ( 6 - 24 ) 
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This can be rewritten in the form 


where we have set 


X?(z) = (Kf * 0 (z) = J Kf(z - Zl )t(zi)&zi , 


I-2T 

K®(t,x)= / Q(s)K £ (t — s,x) ds . 

Jo 


(6.25) 


(6.26) 


It follows that we can set 

(n<f>I)( 2 ) = x?W, (n„!)(2)=x?(2):=(ff°*{)(2). (6.27) 

The second step is to consider products of terms T and ?. To this end, we observe that 
by definition of the canonical model (see (I3.22p l. 

(n £ 0 V)(z) = Xe(z ) 2 = JJ K e (z - Zi)K £ (z - z 2 )£(zi)£(z2)dzidz2 , (6.28) 


which belongs to the second inhomogeneous Wiener chaos. Renormalisation is needed 
because the right-hand side is known to diverge in the limit e —> 0. However, if one 
introduces the Wick product 

£(~i) o £(z 2 ) = £(zi)£(z 2 ) - S(z! - z 2 ) , (6.29) 


then one obtains 


(n £ 0 \r)( z ) = jj K £ (z - zi)K £ (z - Z2)£(zi)o£(z 2 )d Zl dz2 + j K e (z - Zl ) 2 dzi . (6.30) 

The first term on the right-hand side belongs to the second homogeneous Wiener chaos, 
with 

(C W^ 2 W){z)){z 1 ,Z 2 ) = K £ (z-z 1 )K £ (z-z 2 ) , (6.31) 

and is known to converge as e —>• 0. We can thus define 


(n|) c) V )(z) =11 K £ (z - Zl )K e {z - z 2 )i{zi) o£(z 2 ) dzi dz 2 =: xT( z ) > 


(6.32) 


which is indeed compatible with ()6.13bll . provided we set 

C'i(e) = J K e (z - zi ) 2 dz\ = j K e (—zi) 2 dzi . 


(6.33) 


The limiting model (no V )(z) is then naturally defined by the same expression as in (|6.32l) . 
but with K e replaced by K. 

In order to deal with the term 'V ', the relevant Wick product formula is 


f(zi) o£(z 2 ) o£(z 3 ) = £(zi)£(z 2 )£(z 3 ) 

- i{zi)5{z 2 - z 3 ) - £(z 2 )S(z 3 - zi) - £(z 3 )5(r - z 2 ) • (6.34) 

In view of (I6.13cl) . this is indeed compatible with the definition 

(n^*^)^) = JJJ K e{z-zi)K e {z- z 2 )K £ (z-z 3 )^(zi)o^(z2)o^(z 3 ) d Zl dz 2 dz 3 , (6.35) 
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without the need to introduce a further renormalisation constant. 

Before being able to deal with terms such as *Y*, we have to compute the model 
associated with V. Applying (13.271) (in which the sum only contains the term t = 0) and 
the last relation in (|3.28l) . we obtain 

(n £ 0 Y )(z) = f [k(z-zi)-k(- Z i )\ (ngvx^Odsi 
= J [K(z - z x )~ K(—zi)\ xf(z i) dzi 

(6.36) 

where we have used the fact that the integrals of K(z — •) and K are equal to cancel out 
the term C\ . It thus follows from (I6.13cp that 

(nj° ^)(s) = [{K * xf)(z) - (K * xf m]xf (z) -C 2 (&). (6.37) 

Similar expressions hold for the other terms of J- 2 . In order to compute the Wiener chaos 
decomposition of expressions such as (16.371) . we will need the identity 

(£(zi) o £(z 2 )) (£(z 3 ) o£(z 4 )) = £(zi) o £{z 2 ) o f (z 3 ) o U 

+ £(*i) o^{z 3 )5{z 2 - Z A ) + €(zi) o £(z4)8{z 2 - z 3 ) 

+ {,(z 2 ) o^(z 3 )S(zi - z 4 ) +(,(z 2 ) 0 ^ 4 ) 5 ^ - z 3 ) 

+ S(zi - z 3 )5{z 2 - z A ) + 8(zi - z 4 )5(z 2 - z 3 ) , (6.38) 


which follows from m Lemma 10.3] and the Wick rule // ; (/) 0 /^( 9 ) = h+e(f ®g)- It turns 
out that the important term regarding renormalisation is the contribution to the zeroth 
Wiener chaos, resulting from the last line in (16.381) . which is given in the case r = *Y* by 


(W( £; °)-Y-)( 2 ) = 2 


+ 2 


- C 2 (*&) . 


K(z - z)K £ (z - Z\)K e (z - z 2 )K £ (z - Zi)K e (z - z 2 ) dzdz 3 dz 2 

K(-z)K e (z - zi)K e {z - z 2 )K £ (z - zi)K e (z - z 2 ) dzdz 3 dz 2 

(6.39) 


As shown in m Thm. 10.22], only the first integral on the right-hand side (which is 
independent of z by translation invariance) diverges as e —>• 0 , which imposes a choice 
for C 2 (&). In order to represent this quantity and related ones appearing when dealing 
with the other terms, we introduce the notations 


In particular, we have 


Qo( z ) = J K £ (zi)K £ (z 1 - z) dzi , 


Ql(z) = j K £ ( Zl )K?{z 1 - Z) dZ! , 


QH Z ) = f K?(z 1 )K?(z 1 -z)dz 1 . 

(6.40) 

C 2 {&) = 2 J K(z)Q £ 0 {z) 2 dz . 

(6.41) 
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Before proceeding to the convergence result, we collect a number of bounds on integrals 
involving the kernels K e and I\9. In what follows we will sometimes write Kq in place 
of K as well as K® in place of K®, while |x| = \x\\ + \x 2 \ + |a^ 3 1 denotes the i 1 -norm of 
x £ M 3 , so that ||(i,x)|| s = \t\ 4 / 2 + |x|, cf (13.161) . Note however that \x\ and the Euclidean 
norm ||x|| are equivalent. 

Lemma 6.2. The following bounds hold for all z = (t,x) £ K 4 and all e £ [0,1]: 

'^>1 * ^ * WTe ■ (642) 

and 

I K £ (z ) 2 dz < ^ , J K £ (z)k9(z) dz<l, j Kf{zf dz<l. (6.43) 
Furthermore, one has 

leswi < ii- 4 —, ic 5 wi<i, iciwi < i. ( 6 . 44 ) 

Finally, for i,j £ {0,1, 2} let 

hj (e) = J K (z ) Qf (z) Qj (z) dz . (6.45) 

Then Ioo(e) x |loge|, while all other Iij{e ) are bounded uniformly in e. 

The proof is given in Appendix [Bj It is important to emphasize that Lemma 16.21 is 
the key result to show that certain terms do not have to be renormalised. For example, 
consider the difference between the first bound in (16.431) and the second and third bounds. 
This explains why certain symbols involving the ^-operator do not have to be renormalised 
while the same symbol without an additional ^-operator has to be renormalised; see also 
Appendix [B] for a detailed computation. 

Lemma 6.3. For any e, 8 £ (0,1] and any (t,x) £ M 4 , one has the bound 

\KQ(t,x)-K?(t,x)\<-~L_. (6.46) 

As a consequence, 

j (K?(z - Zl ) - kQ{z - Zl )) {Kf(z - Zl ) - K^(z - Zl )) dz! < e 2e (6.47) 

holds for all 9 £ (0, |) and all z, z £ M 4 . 

We give the proof in Appendix O 

We can now state the main result of this section, which is an adaptation of m 
Thru. 10.22], 

Proposition 6.4 (Convergence of the renormalised models). There exists a random model 
Z, independent of the choice of mollifier g, and a family of renormalisation transforma¬ 
tions M £ £ 91 such that for any 9 < — | — ao = k (recall from (13.61) the definition of the 
noise regularity ao ), any compact set and any 7 < £, one has 

E||| M £ Z £ ; Z|| 7i s < (6.48) 

uniformly over e £ (0,1]. There Q is such that (13.25[) holds for all polynomials P of parabolic 
degree less or equal f. 
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Proof: As already stated, the proof follows along the lines of [19, Thm. 10.22], with some 
changes due to the presence of the additional integration operator £. Thus we will only 
comment on those parts of the proof that need to be adapted. 

The result will follow if we are able to choose the constants C\ (r) and C 2 [r) in (16.111) 
for each e G (0,1] in such a way that the bounds (16.161) hold for all r G equivalently, 
the bounds (16.201) should hold for all these r. The limiting model Z will be constructed 
step by step during this process. 

The first step is to set 

(n 0 3)(z) = £(*), 

and the fact that the required bounds are satisfied is a consequence of Proposition 9.5 
in m (see also Lemma 10.2). 

In the case of r = T, the limiting model is defined by (16.231) , and the conclusion follows 
as shown in [ 191 Thm. 10.22] by applying the bounds on convolutions in [T9’, Lemmas 10.14 
and 10.17]. We thus turn to the case r = ?. By (16.271) . we have 

((^+(z))(z 1 ) = i+(z-zi), 

((n ,(1) ?)(z))(zi) = K q (z - Zi) . 


The first bound in (16.201) then follows from the fact that 

((W ( 1 ) ?)(z),(W ( 1 ) ?)(z)) = J K q (z- Z l )K Q {z-z l )dz 1 = Q° 2 {z-z) 


(via the change of variables z\ H > z — z±). The right-hand side being bounded as shown in 
Lemma [6721 the bound (|6.20[) is satisfied for A = 0 provided R ^ —2| ? | s = 1 + 2k. In order 
to prove the second bound in (16.201) . we use the bound (16.471) in Lemma 16.31 which yields 


((<5W (e;1) ?)(f, x), (dW (e;1) ?)(£, x)) < e 29 

for any 6 G (0, ~). Next we consider the term V. Then a decomposition analogous to the 
one given in (16.30[) shows that we have 

((WWV)(Z))( Z1 ,Z 2 ) = K e (z-z 1 )K?(z-z 2 ) , 

(VV^ 0) V)(z) = J Kei-zJKQi-zJdzx-Ciie) , 


where we have used the change of variables z\ i— > z\ — z. Note that the second line is 
independent of z, as a consequence of translation invariance, so that we can drop the 
argument z from the notation. The integral in the second line is equal to Q\ (0), which 
is bounded uniformly in e by Lemma 16.21 We may thus simply choose C\ (e) = 0, which 
amounts to setting 


(nj° v)(*) = Jl K £ (z - Z!)K?(z - z 2 )Z(z l ) o t(z 2 ) dzi dz 2 + V 
= Xe(z)xf(z) ■ 


Indeed, the first bound in (16.201) for k = 0 is satisfied with A = 0, provided R ^ —2| V | s = 
2+4 k. Note that being allowed to set C((e) = 0 here is precisely the point where we use the 
fact that certain terms do not have to be renormalised due to the bounds in Lemma [ 6 .21 Of 
course, another possibility would be to set C[(e) = V, which would result in IIq V 
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belonging to the second homogeneous Wiener chaos. Since Qf(0) is bounded, however, 
the choice of C[ (e) does not really matter, since it only results in a shift of the random 
variable defining Z. Whatever the choice, the bounds (16.201) can be checked by similar 
arguments as in the previous case. 

The symbols V, ^^ as well as V Xj and VX, can be treated in pretty much 
the same way. It thus remains to consider the symbols in the families of 3*, 3 and 3*. 
As shown in the proof of [191 Thm. 10.22], the symbol 3* already fulfils the bounds (16.201) 
as a consequence of the choice of C\{e). Since the kernel K® is less singular than K, the 
same conclusion holds for all symbols of the same family. 

The symbol 3 > requires the introduction of a further renornralisation constant 3 > ) 
which has order |loge|. We thus have to check what happens to the other symbols in the 
same family. Consider for instance the case of *3*. A computation similar to the one 
made in (16.361) and (|6.37l) yields the expression 

(ng*3*)(z) = [(K * (xe o - (K * (xe<>xf))(0)\xf(z) - C 2 (*3*) . 

Writing out the integrals and applying (|6.38|) . we can decompose this expression into a 
sum of terms in the Wiener chaoses of order 0, 2 and 4. In particular, the term in the 0th 
Wiener chaos is given by 

(W( e;0 ) *&)(z) = 2 jjj K(z — z)K £ (z — z\)K® (z — z 2 )A' e (z — z\)K e (z — 2 2 ) dzdz\ dz2 
+ 2 JJJ K(-z)K e (z - zi)K®(z - z 2 )K £ (z - z 1 )K e (z - z 2 )dzdz 1 dz 2 

- c 2 (3*) • 


The first integral on the right-hand side does not depend on z, as can be seen by the 
change of variables (z,z\,Z 2 ) (z + z, z± + z,Z 2 + z). As a consequence, it is equal to 
2/oi(e), and is thus uniformly bounded in £ as shown in Lemma 16.21 The second integral 
is also bounded, so that we may choose C2(*3 > ) = 0 (or any other finite value). The 
terms in the other two Wiener chaoses are already bounded in the case of *3*, so they 
remain bounded in the present case. It is now quite obvious that the other symbols of the 
family do not need to be renornralised either, as the relevant integrals are obtained from 
the above ones by substituting the appropriate kernels K by K®. 

Finally, we have to deal with symbols in the family of 3*. However, it is shown in the 
proof of [19, Thm. 10.22] that the constant C , 2(*3*) suffices to renormalise 3 * as well, so 
the situation is completely analogous for the other symbols of the family. □ 

6.3 Computation of the renormalised equations 

In order to complete the proof of the main results, it remains to derive the classical equation 
with mollified noise that corresponds to the renormalised models. In other words, for a 
given M = M £ £ 94, we have to characterise the solution map Sfa satisfying 

TI m Sm{u 0 ,vq,MZ £ ) =Sm(u 0 ,v 0 ,^ £ ) , (6.49) 

where 7 Z M is the reconstruction operator of the model (= n( £ )) and Sjfa is the map 
giving the fixed point of the renormalised equation in terms of initial data and forcing 
term. Relation ()6.49|) can also be represented as the following commutative diagram: 
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(u 0 ,v 0l MZ £ ) 

MT 


( Um>V*m) 



n M 

{uo,v 0 ,£ £ ) I-—-> ( u £ ,v £ ) 

S L m 

In order to determine the classical equation, we will need the following algebraic result. 
Lemma 6.5. For all r € F o, one has 

(n fr){z) = (n z Mt)(z) . (6.50) 

Proof: The claim can be checked directly by comparing the expressions given in (16.121) 
and (16.131) . More generally, the first relation in (16.41) shows that (16.50[) holds for all t 
such that A AI (r) = Mr < 8 > 1, which happens to be the case for all symbols excepts those 
of the families and *$*. In the exceptional cases, the relation has to be checked 

by an explicit computation; in general it holds only when the model is evaluated at the 
base point, i.e. we may have {H¥t){z) (II z Mt)(z) if z z. For instance, in the case 
r = , the action of the renormalisation map is given by 

= % - 2c[ i* - c'{V . 

One can then check that in the case r = ^, ()6.5D is satisfied by setting A M ^ = M ^ < 8 > 1 
and MJ (^) = J (M ^). Using this and the definition of the map A, a somewhat lengthy 
computation shows that 

satisfies the second relation in (16.51) when r = ^*. By (I6.4j) , this implies 

nf t- = (n 2 ®/ 2 )A M t- =n 2 Mt* +2C[n z (lx i )(f z , jj) + c'fn z (tx i ){f z ,j i l ), 

where we used the fact that (/ 2 ,1) = 1. Since (II z $Xi)(z) = (n 2 ?)(z)(ak — Xi ) by ()3.22j) . 
we see that the last two terms on the right-hand side indeed vanish when z = z. □ 

As a consequence of (I6.50p . when / is a function, as is the case for e > 0, we have 

(H M f)(z) = (Bf/(*))(*) = ( U z Mf(z))(z) = (KMf)(z) . (6.51) 

This directly implies the following result, which is the equivalent of Proposition 15.11 for 
the renormalised equation. 

Proposition 6.6 (Renormalised equation). Assume that e > 0 and that (15.11) admits a 
fixed point (U,V) where U, V £ T> 7,r> for some 7,7 G M. Then ( ii,v) = (7Z M U,7Z M V) 
satisfies the classical SPDE 

d t u = Au + F(u, v) + ffi , 

dpu = Aiii + A 2 V , (6.52) 

provided F is such that 

F(MU , MV) = MF(U, V ) + g(U, V ) , (6.53) 

with g having only components of strictly positive homogeneity. 
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Proof: Applying M to the expansions of U and V given in (15.61) . we see that MU differs 
from U and MV differs from V only by components of homogeneity | —n at least. Applying 
1Z to MU and MV and using (16.511) . we obtain that as in the proof of Proposition 15.11 


u{z) = Xe{z) + V(z) , 
v(z) = x?(z) + fi(z) . 

On the other hand, applying M to the right-hand side of (15.31) and using the fact that M 
commutes with X and £ and leaves polynomials invariant, we get 

MU = 1R+ [H + MF(U, V)] + R(MU, MV) + Gu 0 , 

MV = SR+MU _ + QR+MU+ + Qv 0 . 

Applying the reconstruction operator, we thus obtain 

u(t, x) = G * R + [£ £ (f, x) + ( 1ZMF(U , V))(t, x)] + Guo(t, x ) , 

v(t,x)= / Q(t — s)v(s, x) ds + Qvo(t, x) . 

Jo 

The claim then follows from the fact that (|6.53l) implies 

(1ZMF(U, V))(z) = ( KF{MU,MV))(z ) = F(u(z),v{z)) , 

where the last equality follows from the same argument as the one used in the proof of 
Proposition 15.11 □ 

Finally, we compute the expression of F. To simplify the computations, we assume that 
Cl = C'i(V) and C 2 = C , 2 (*v^') are the only non-zero renormalisation constants, which 
is a possible choice as seen in the proof of Proposition 16.41 This is also the place where 
the assumption 72 = 0 turns out to be necessary in order to obtain a simple renormalised 
equation. 

Proposition 6.7 (Expression for the renormalised F). Let F be the cubic polynomial 
given in ([53D with either d = 2, or d = 3 and 72 = 0, and let M be the renormalisation 
transformation defined in (16.111) , where the only nonzero constants are C\ = C\ (V) and 
C 2 = C 2 { : &). Then 

F(u , v) = F(u, v ) — co(e) — ci(e)it — c 2 (s)v , (6.54) 

where the coefficients are given by 

c o(e) = fii{Ci + 371 C 2 ) , 
ci(e) = 371 (Ci + 371 C 2 ) , 

02 (e) = 72^1 • (6.55) 

Proof: The fact that (U, V) differs from (MU, MV) only by terms of order | — k implies 

F(U, V) = MF(U, V) + e'(U, V) , 

where all components of g'(U, V) have strictly positive homogeneity. Thus we only have to 
compute the difference MF(U,V) — F(U,V), which reduces to determining the terms of 
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F(U,V ) which are modified by the renormalisation map, or, more precisely, which terms 
of homogeneity less or equal 0 appear when applying M. 

The computation is straightforward, though somewhat lengthy. In practice, it suffices 
to expand the quadratic terms of F(U, V) up to order \ — k, and cubic terms up to order 
0 — 2k. To explain where the condition 72 = 0 arises if d = 3, we allow for general 72 € M 
for the time being. 

Consider first the quadratic part F 2 (U. V) = I3\U 2 + /3 2 UV + P 3 V 2 of F(U, V). It turns 
out that when the renormalisation map M is applied, only the term proportional to V 
of F 2 (U,V) yields a contribution of not strictly positive homogeneity. Namely, one has 

MF 2 (U, V) - F 2 (U, V ) = -PxCxl + 02(17, V) , 

where q 2 ( U. V ) has strictly positive homogeneity. It remains to determine the contribution 
of the cubic part F 3 (U. V) = 71 U 3 + 72 U 2 V + 'y 3 UV 2 + 74 V 3 . Substituting U and V by 
their expansion (|5.6j) and sorting terms by homogeneity yields an expression of the form 

F 3 (U, V ) = F 3 ._s_(U, V) + F 3 ._i(U, V ) + F 3 ._!(17, V) + F 3 . 0 (U, V ) + q 3 (U, V ) , 

where each F 3 .g contains only terms of homogeneity /3 — 0{k) and the remainder g 3 has 
strictly positive homogeneity. We treat all of these terms separately. For instance, we 
have F 3 ._3(U, V) = 71 + 72 ^ + 73 *^ + 74 ^, where the first two terms have to be 

renormalised, yielding 

MF 3 ._3{U,V) - F 3 ._3(U,V) = - 371 C 1 T - 72 C 1 ? + e 3 ._3(u,v). 

The term F 3 --i(U,V) yields a contribution 

MF 3 .-i(U, V) - F 3 ._i(U. V ) = -3 7 iCi^1 - l2 C^l + 0 3 ; _i(Z7, V) 

stemming from the terms V. Consider next the term F 3 .q(U. V). It contains 2 terms 
proportional to , whose renormalisation yields the contribution 

MF 3 ,o(U, V) - F 3 . 0 (U, V) = - [37161 + 7261] C 2 1 + g 3;0 (U, V ) . 

Finally, in the term F 3 ._i(U,V), it is the terms in •s and that need to be renor¬ 
malised. One obtains 

MF 3 ._i(U, V ) - F 3 ._i(U, V) = - [971 ai + 37 2 ai] C 2 T 

— [ 971^2 + 372 ^ 2 ] C 2 ? + g 3 ._i (U,V) . 

In the case 72 = 0, using ai = 7 i and replacing the 6 * by their expressions (15.71) . and adding 
the last four expressions, we see that we can always factor out the expression ^> 1 + T, which 
coincides with U up to terms of strictly positive homogeneity. The result follows. 

In case 72 / 0, there are additional terms such as 72 C 261 I, which do not admit a 
simple expression in terms of U and V, so that the renormalised equation has no simple 
closed-form expression. However when d = 2, then we can choose C 2 = 0 because the 
symbol has strictly positive homogeneity and does not have to be renormalised. Then 
we can also factor out the quantity ipl + ?, yielding the expression for C 2 (e) in (16.551) . □ 
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7 Proof of the main results 


Since Theorem 12.11 is a particular case of Theorem 12.21 we proceed directly to proving the 
latter result. 

Proof of Theorem 111. HI Consider first the situation on a fixed, sufficiently small time in¬ 
terval [0,T], Fix e > 0 and let 

(U £ ,V £ )=S T (u 0 ,v 0 ,M e Z £ ) , 

(U,V) = S T {u 0 ,v 0 ,Z) , 

where St is the solution map given in Proposition [53] (we have suppressed the dependence 
on W, which is the same for both fixed points). For any 5 > 0, the bound (15.161) . Markov’s 
inequality and Proposition 16.41 imply 

Pjlll U e ; U\\ 1iV] t + \\V £ ; V\\^, T > 0)5'} < ^||| M e Z £ ; Z||| 27+Q;0 < j . 

It follows that (U £ , V £ ) converges to (U, V) in probability as e —> 0. 

Next, let (u,v) = (1ZU,1ZV) and (u £ ,v £ ) = (1ZU £ ,1ZV £ ) (where we do not indicate the 
model-dependence of the reconstruction operators 1Z). The fact that (u £ ,v £ ) converges to 
(u, v ) in probability as e —> 0 is a consequence of the bound [T9], (3.4)] in the reconstruction 
theorem, combined with the definitions (13.201) of ||II|| 7; ^ and (I4.28|) of ||| • ||| 7)?7; £. 

The convergence result can be extended from a fixed time interval [0, T] to the random 
time interval up to the first exit from a ball of radius L in exactly the same way as in m 
Section 7.3]. In particular, the required continuous dependence on the model is proved 
in PS Cor. 7.12], 

Proposition 16.61 shows that (u £ ,v £ ) satisfies the system (12.111) . where the explicit ex¬ 
pression of F is a consequence of Proposition [677] The expressions (I2.5|) and (12.121) for the 
renormalisation constants in dimension d = 3 follow from (I6.55[) . where the only nonzero 
terms are C\ (s). defined in (16.331) . and C 2 (£), defined in (j6.41|) . 

In the case d = 2, Table |T| shows that only the terms T, V and ^7 need to be 
renormalised, to that the expression m for the renormalisation constant follows by 
taking (^(e) = 0 in (16.551) . □ 

Proof of Theorem 1 2.ZA The only difference with the previous case is that we have to take 
into account the presence of the n different operators Si, i = 1 ,...,n, introduced in 
Remark 14.21 This means that the fixed-point equation for V in (15.11) is replaced by n 
different equations of the form 

Vi = IC^R+U + QiV o , i = 1 ,..., n . 

Since n is finite, it is straightforward to check that the proofs of Proposition 15.11 and 
Proposition 15.51 carry over to this situation. 

It remains to check the expressions (|2.18p for the renormalisation constants. To com¬ 
pute the diverging part of F(U,V) — F(U,V), it is in fact sufficient to determine the 
coefficients of V, and •3* in F(JJ,V ), since these are the only ones that need to be 
renormalised. These can only result from the terms j3\U 2 + 71 U 3 in F(U,V), since all 
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other terms will contain at least one of the £j . It is now sufficient to check that U has the 
same expansion as in (1 5.6 [) . except that the first relation in (15.71) becomes 


bi = /3i + 3(^71 + ^2 ^72,1 , 

i=l 

where ipi is the coefficient of 1 in V)- The result then follows as in the proof of Proposi¬ 
tion [6/71 taking into account the modified value of b\. □ 

A Proof of Lemma 14.8 

The proof is by induction on the size of the regularity structure, which is constructed 
recursively, starting with S, 1 and the 7Q, and applying the operators I and £ as well 
as multiplication to already existing symbols. We show that the properties (|4.21|) hold 
on the polynomial regularity structure and for r = E, and that they are stable under the 
algebraic operations extending this structure. 

1. The claim certainly holds on the polynomial part of the regularity structure, since 

0 ni +h x k )(z + h) = {z + h-z-h) k = (n %x k ){z ), 
n +h ,- z+h X k = (X-z-h + z + h) k = T\ z X k , 

and 

r l +h - z+h x k = (id®^ z+hrz+h )Ax k = x k + itf z+h - z+h ,x k ), 

n,z Xk = (Id®7,%-)AX fc = + l^ £ z - zl X k ) . 

2. Using the fact that A(H) = E ® 1, it is also immediate to check that (|4.21|) holds for 
r = E (this reflects translation invariance of the noise). 

3. We show that if the relations (14.211) hold for r, then they also hold with r and f 
replaced by It and J^t. Indeed, 

{f £ z+h , Jkr) = - J D k K(z + h- z)(U £ z+h r){z) dz 
D k K(z-z)(U £ z+h r)(~z + h)dz 

D k K{z- z)(U e z r)(z)dz 
= (. fz,JkT) , 

and 

(n e z+h X T ){z + h)= f K(z + h — z')(n £ z+h T)(z') dz' + p-(/!+*> Jt T ) 

J l 

K(z - z){I\. £ z t){z) dz + J2 ^ 

l 

= (n:-xr)(z). 
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Furthermore, it is shown in the proof of [191 Thm. 8.24] that 

r !+M+fc( Xr ) = (M®7 t+h,z+h)P ® Id ) Ar + ^ ^ ('fz+hj+h, 3t r ) > 

where E M d+1 has components (z 7 )j = ~{ r Yz+h z+h^i)- The first term on the 
right-hand side is equal to XT% +h z+h (r). The induction hypothesis implies that all 
terms are equal to their value for h = 0, proving the required identity for T £ z2 (It). 
The identity for 7 zzi^k T ) follows in a similar manner. 

4. In a similar way, using the definition (14.61) of TT z (£t). it is straightforward to check 
that if r satisfies (14.211) . then St satisfies (14.211) as well. 

5. If r,a satisfy ([4.211) . then 

(U s z+h Ta){z + h) = (n % +h r)(z + h)(U £ z+h a)(z + h) = (TT z t){z)(TT z (t){z) = (II £ z t<t)(z) . 
Furthermore, 

r *+M+ft( r<r ) = ( Id ®7 '* +m+ JA(tct) = ( Id ®7 £ zz) a ( t °) = r ^(™) • 

Here we have used the fact that the induction hypothesis applies because the second 
component of A(r) always depends only on already computed quantities, cf. Table 1. 
The translation invariance of 7fj(' rcr ) follows. □ 


B Proof of Lemma 16.2 


It follows from the definition of K that it satisfies \K(z)\ < Hz]!” 3 . Computing the 
derivatives of K , one sees that it is singular of order —3 in the sense of [T9l Def. 10 . 12 ], 
Thus the first bound in (16.421) follows from [_19i Lemma 10.17]. 

The second bound in ([6.421) is a crucial estimate, which allows us to avoid renormalisa- 
tion of certain symbols involving the operator S. Therefore, we provide additional details 
in the proof of this bound. Recalling the definitions of K® and Q we have 

r r2T/\{t+e 2 ) 

K®(t,x)= / Q(s, x)K e (t — s, x) ds = / Q(s, x)K e {t — s, x) ds 

Jr Jo 

since K e (t,x ) is supported on {f > — e 2 }. By a change of variable s i->- t — s in the last 
integral and using the boundedness of Q on [0, T] for a final time T > 0, it follows that 

K Q (t, x) x f ^ K e (s,x) ds, 


where we recall that the notation / x g indicates that we have both /<§ and g < /. 
Therefore, combining the last characterisation of K®(t,x) with the first bound in ([6.421) 
yields 


\K?(t,x)\ < 
< 


L 


i 


_ e 2 (| S |V 2 + |x| + e ) 3 

KM+e ) 2 l 

£ 2 (M+e ) 3 


ds 


ds + 


L 


(M+e ) 2 


1 

^|372 


ds 


< 


1 

x| + £ 
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and this concludes the proof of the second bound in (16.421) . The first bound in (16.431) 
follows from 


J K '^ dz ~ J m^rw dzi 


1 rl 


r 2 dr 


o Jo (t 1/2 + r + e) 6 


dt 


1 


The second bound in (16.431) is again important, since it leads to terms which involve 6 
but which do not have to be renornralised, so we provide more details in the proof of this 
bound. Note that 


J K £ (t, x)K® (t , x) dx < 


1 


1 


< 

r^j 


< 

rsj 


(\t\ 1 / 2 + |x| + e) 3 (|x| + e) 

1 1 


(|t| 1/2 + (Ml +e) 3 (||x|| +e) 


dx , 
dx 


1 


1 


(|t|V 2 + r + e) 3 (r + e) 


sin^i)?’ 2 d 6 dr 


using the equivalence of norms and spherical coordinates (r, 6) = (r, 61,62) in 
fore, we find the estimate 


There- 


K e {t,x)K®{t,x)dx < 


r 2 dr 


1 


(r + {t] 1 / 2 + e) 3 (r + e) \t\ 1 / 2 + £ 


Integrating over time then yields 


K £ (t,x)K®(t,x) dx df < 


L £ 2 |f|V 2 + e 


dt < 1 


which proves the second bound in ()6.43l) . The third bound in ()fi.44[) follows from the fact 
that 

r 2 dr 


J K ? 4 x) 2 dx < J 


1 


(r + e) 2 

To obtain the bound on Qo, we distinguish the regimes ||z|| s ^ e and ||z|| s > e. If ||z|| s ^ e, 
the same computation as the one made for the integral of K e (z) 2 shows that |Qo(z)| 1$ 1 /e, 
while for ||z|| s > e, we can use m Lemma 10.14] to obtain 


I Go (*) I 


< 


kills Ik-21II! 


d*i< 


Using spherical coordinates, one can show that for all a, b ^ 0 and all n ^ 1, one has 


4 := 


1 


1 


■ dx 1 


r 2 dr 


(|xi| + a) n \x- xi\ + b 
From this we easily get that whenever a ^ b ^ 0, 


/ 0 (r + a) n (r + |x| + b) 


h<l + \x\ + a 2 , I 3 < 


1 + | log a | 
\x\ + a 


and thus 


mt,x)\ 


< 

r^j 


- 1 1 + |log(t 1 / 2 + e)| 


l 0 |x|+U/ 2 + e 
I Qv(t,x)\ < 1 + |x| +e 2 < 1 . 


d* < 1 , 
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It remains to bound the terms /^(e). It is in fact sufficient to check that 


1 rl 


loo(z) 


-foi(e) 


r 2 dr 


/o Jo (t 1 / 2 + r) 3 (t 1 / 2 + r + e) 2 
"* / ' 1 r 2 dr 


di < I log e | , 



o Jo (t 1 / 2 + r) 3 (t 1 / 2 + r + e) 


dt < 1 , 


since the other terms can be bounded above by Elementary computations show 

that these bounds indeed hold true. □ 


C Proof of Lemma 16.3 

Let us recall Definition 10.12 in [19]. Given a scaling s, a smooth function K : M rf+1 \{0} —> 
M is said to be (singular) of order ( if for every sufficiently small multiindex k , one has 

\D k K(z)\ < ||z||s \fz such that ||z|| s ^ 1 . 


More precisely, we require that there exists m ^ 0 such that 

|||^IC;m : = sup sup \\z\\[ k ^ C \D k K(z)\ < oo . 

k: 2gK d +! 


Note that in QE Def 10.12], the second supremum indeed runs over the whole space, but 
in practice we will only apply it to compactly supported functions. Lemma 10.17 in m 
states that if K is of order £ £ (—|s|,0), then K e = K * g £ has bounded derivatives of all 
orders, satisfying 

io^wi <(ii2ii.+n c - w -iA'ii iW ,. 

Furthermore, for all £ £ [£ — 1,£) and m ^ 0, one has 

IK - K £ | c - ;m < e^lKl c . m+max{gi} . (C.l) 

Consider now the quantity 

K®(t, x) — K®(t, x) = f [K(s, x) — K £ (s, x)] Q(t — s) ds 

Jo 

(in fact, the lower bound can be replaced by 0 V (t — 2 T)). We have 

\K Q (t,x) - K$(t,x)\ < [ |iL(s,x) -K e (s,x )\ds . 

Jo 

We already know that K is singular of order —3, and thus |||/i|||_ 3 ;m is bounded for any 
77 i. In particular, we have by (jC.lD with m = 0 the bound 


P'-AL|| c - ;0 < e - 3 -q||Al-3 ; 2 

for any ( £ [—4, —3), which implies 


| K(s,x) - K £ (s,x )| ^ 


\\k-k £ \i^ 0 ^ g -3-c 


Nils C 


Mis 


l*l-3;2 • 
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Taking £ = —3 — 6 for some 9 £ (0,1], we obtain 


\K(s,x) - I< £ (s,x)\ 


< 


ui|3+0 
\^\\ s 


l*1-3;2 


and thus, since |||it ||- 3 ; 2 < 1, 

\K Q (t,x)-K?(t,x)\<s e g(3+e)/2 d ^ |a 


< 


,|3 +0 - | x |l+0 > 


proving (16.4611 . To prove (16.4711 . we use spherical coordinates as in the proof of Lemma [6.21 
to get, for any constant c of order 1, 


dx\ 


sin <f> d <f> 


'|au|<c \x 1 \l+°\x-X 1 \l+<> 


r 1 r 2 r 

Jo rl+e Jo (x 2 + M 2 — 2r|x| cos 0)( 1+0 )/ 2 

f 1 (r + |x|) 1-0 — |r — |x Ml 9 


dr 


< 


L 


2|x| \ x \i-e 


2\x\r e 
dr 


dr 


+ 


f 

J 2\x\ 


|x|r" 
i r i-e 


ir 


l-e 


1 + — - 1 - — 


i-0' 


dr 




'2|a:| 


dr 

720 


<1 


for \x\ < 1, provided 29 < 1. The result follows since t >->• K®(t,x) has compact support. 


□ 
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